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ABSTRACT 


This it th« flntl report submitted by SAO nader contract 
NAS8-3S03<$« "Investigation of Electrodynamlo Stabilization and Control 
of Long Orbiting Tethers," and covers the period 13 September 1982 
through 12 September 1983. It details the results of research 
undertaken by SAO on three topics under the general area of the 
dynamics of spacebome tethers: 

1) SAO reviewed existing dynamics models and 
software and made recosusendatlons for 
future software development. 

2) SAO developed a computationally efficient 
high-resolution tether model and performed 
studies of the dynamic behavior of a 
broken tether. 

3) SAO studied three advanced applications of the tether: 

a) Oamping of longitudinal tether oscillations by 
means of a reel motor control algorithm and 
use of the technique In a payload orbital 
transfer problem, 

b) Payload a cquisition by a Sjjtce Station, 
c^ Orbital p'ump in g . 


The authors of this report are Or. Giuseppe Colombo, Principal 
Investigator: Mr. Oavid A. Arnold and Or. Gordon H. Gullahorn, 

Analysts: and Mr. Richard S. Taylor. Program Manager. Mr. Arnold was 

also a Co-Investigator. 

Or. Colombo is a Visiting Scientist at SAO from the University of 
Padova, Italy. 
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T«ther Dynamics Software Applications Matrix 
Tether Dynamics Software Coapnter Facilities and 

Langnages in nse. 

Figure 4-1. A SITYHOO^ run showing tether recoil after a break at 
200 meters of a 100-km tether of 0.2-cm Kevlar with a 300-kg 
snbsatelllte. The figure is in the familiar side-view form. The 
in-plane component is forward (in the direction of Shuttle motion) to 
the left. Configurations are plotted at 2-second intervals, as they 
are in Figure 4-2 to Figure 4-5. 

Figure 4-2. SLACK2 was used to run a close analog of the case in 
Figure 4-1. The tether was divided into four segmeuts. Note the boom 
(solid liue) deployed straight up, with the boom tip as "laass 6 ." The 
curvature in some dashed lines is an artifact of the plotting process. 
The vertical scale is not quite the same as Figure 4-1. In computing 
the recoil velocity, the tether mass was set to zero (see text). 

Figure 4-3. As in Figure 4-2, with 40 tether segments. For 
large numbers of segments, the plotting program was modified to omit 
plotting the "mass points" other than the boom tip and Shuttle. 

Figure 4-4. In plane vs. radial behavior of the tether plotted 
every two seconds for a total period of 240 seconds after a break for 
the situation modelled in Figure 4-2 and Figure 4-3 except that tether 
mass is included in computing the initial recoil velocity and the boom 
is deployed 30*^ forward of the local vertical. The first 52 seconds 
of the run is shown above (a); the entire run is shown on the next 
page (b) . Note the tether oscillations induced by the motion of the 
boom and the overall behavior similar to damped oscillations. 

Figure 4-5. A series of runs for the same case as in Figure 4-4, 
with 10-, 20-, 30-, 40- and 5o-tsther segments. All runs are for 38 
seconds with output at 1-second Intervals. Note the qualitatively 
different results obtained with 10 and 20 segments in (a) and (b) from 
those with 30, 40 and 50 segments in (c) through (e). 

Figure 4-d. Timing results for the four runs shown in Figure 
4-5, plotted as cumulative results for the full 38 seconds of each 
run: (a) The ratio of the computation (CPU) time to the physical time 

modelled; log-log plot vs. the number of segments: and. (b) The 

ratio of the number of bounces to the physical time: log-log plot vs. 

the number of segments. Also plotted in (b) is the CPF time taken to 
compute the average bounce, which is seen to increase much more slowly 
with increasing resolution. 
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Figure 3-1. 
Figure 3-2 . 


Figare 4-7. The final configurations of three separate runs are 
plotted side-by-side. The physical sitnation modeled in each ran is 
the same (as in Figure 4-4). The tether is divided into 30 segments 
in each case. Each ran is for 40 seconds. The three cases differ 
only in that after an initial partition of the tether was computed, 
the individual segment lengths were varied by small random amounts in 
the range - 1 % to 2%. The lengths were then normalized to give the 
same total length in all oases: (a) shows the configurations after 40 

seconds; and (b) shows the configurations after the full $0 seconds. 
See the remarks in the text for cautions on the interpretation of 
detail in the figures. 

Figure 4-8. Part (a) shows the evolution of one of the cases 
nsed in Figure 4-7, while (b) shows a similar case in which the tether 
was more evenly segmented resulting in less resolution near the boom 
tip. Note the distinctly different behavior of the two cases. In 
(a), a travelling wave was generated, while (b) exhibited no such 
wave. This same distinction is seen between Figures 4-S(a) and (b) 
and Figures 4-5 (c), (d), and (e). The finer resolution results are 
more consistent among themselves and probably closer t" the real 
behavior of the tether. 

Figure 5-1. Behavior of the roots of the characteristic equation 
for fixed values cf T. Note that we plot the negative of the real 
part of each root (£sirs of complex roots are indicated by a dashed 
line) and that for K, 3 > *) nil roots have negative real parts (i.e. 
the system is stable).^ The quantity of greatest Interest is the lower 
line at each value of 8. This root governs the solution which decays 
least rapidly, and we want to choose parameters^K, 8 to optimize this 
worst solution. In choosing actual parameters K, 8 we may not achieve 
the desired dimensionless parameters, so there is the additional 
constraint of robustness: slight errors in the parameters should not 

lead to large decreases in performance. 

Figure 5-2. Tension vs. time for a release study (100-ton 
Shuttle, 9,5-ton payload, 0.5-ton teleoperator, and an 80-km length of 
0.3-cm Sevlar tether) in which the tether is reeled out first to 
relieve tension. The maneuver is given by (deployed length) - 80 km 
A[l-cos(2nt/P) 1 . with amplitude A « -650 meters, period P » 209 sec. 
After 1/2 period, the deployed length is held constant. The payload 
is released at 85 seconds. 

Figure 5-3. Tension vs. time for a release study in which all 
parameters are the same as in Figure 5-2 except that at t > 100 sec 
the sinusoidal maneuver is terminated (held constant) and t^e active 
damping control law_^s initiated (parametersT = 1.8 x 10 dynes; T 
■ 1.4 X 10 cm-sec dynes ; 8 * 0.041 sec '?. 

Figure 5-4, Tension oscillations after payload release, with no 
damping. A reel-in/reel-out sine-law maneuver was used prior to 
release with amplitude 174 meters, period 114 seconds, phase offset 
-n/2. The manenver terminated at 114 seconds and release was at 117 
seconds. We plot only from t * 105 seconds. System: 80-km tether of 

0.3-cm Kevlar; 100-metric-ton Shuttle; 9.5-ton payload; 0.5-ton 
remaining teleoperator. 



Figure 5-5. Tension oscillations after release^ with active 
damping* The system and pre-release maneuver are the same as for 
Figure 5-4, and the damper startedi^at release, t » 117 seconds^ 

1)am£|r parameters: T « 2*5^ a 10 x 10 dynes, K » 1*6 z 10 cm 

sec dyne , p » O.oSl sec 

Figure 5-6. The same as Figure 5-5, except that the damper 
parameter VC and 9 are 75% greater* The response is not significantly 
altered. 

Figure 5-7* The same as Figure 5-5, except that the damper's 
nominal tension T is $3% of that used in Figure 5-5. This figure and 
Figure 5-6 above illustrate the robustness of the damper algorithm; 
its success does not depend on choosing optimized parameters* 

Figure 5-8, Tension variation with a pure reel-out pre-release 
maneuver and active damping after release. Note that the initial 
tension increase is avoided* The system configuration is as in Figure 
5-4. The pre-release sine-law maneuver had an amplitude of 406 
meters, a period of 200 seconds, and a phase offset of -n/2* Release 
is at 70 seconds, and the maneuver terminates at 100 seconds (P^2). 
Active damping begin| at release with^parameters T^ 
dynes, K » 1*6 x 10 cm sec dynes , ^ = 0.041 sec 

Figure 5-9, Tether deployed length as a function of time. 
Components due to sinusoid reel-out maneuver and active damping 
response shown separately. 

Figure 5-10. Reel velocity as a function of time* Components 
due to the sinusoid reel-out maneuver and active damping response 
shown separately. 

Figure 5-11, Tether tension from the initially deployed 
configuration (t ^ 0) through the release and damping (t »1441) until 
t * 2000 seconds* The detailed mission profile is described in the 
text* 


Figure 5-12. For the case in Figure 5-11, the side-view 
(in-plane vs* radial components) of the system is plotted at 10-second 
intervals. Each view is arbitrarily displaced in the in-plane 
direction to avoid overlap, the earliest configurations being at the 
left* The components are relative to the Shuttle, and uniformly 
spaced constant vectors have been subtracted from each mass's position 
(effectively "cutting out" a portion of tether) to expand the 
resolution in the radial direction. 

Figure 5-13* The tether tension is plotted near the time of 
release, with an expanded time scale. The pre-release maneuver and 
active damping are more readily seen. 
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Figure 5-14. A sohematlo view of the various orbits discussed in 
the text. The heavily drawn portions represent actual trajectories. 

Figure S-IS. X vs. ft. X is computed from the equation 
[l-Xl’tl-(X+ft) /ll-al-ft for the parameter ft with range 0 to 1. X is 
the scaled tether length (L/R) and ft is the similarly scaled 
difference in height between the station orbit and the payload parking 
orbit. The ratio of the actual solution X to the first-order 
approximation ft/7 is plotted. 

Figure 5-16. Tether tension vs. time for a case in which the 
payload is acquired with no other maneuver or dr ping. The system 
(1000-metric-ton Space Station, 0.5-ton teleoperator. 44 km of 0.3-cm 
diameter Revlar tether) was initially in equilibrium in a 500-km 
circular orbit, with the tether deployed downward. A 9.5-ton payload 
was acquired with zero relative velocity (payload delivery orbit, 200 
to 456 km) at t > 10 sec. Note the oscillation with period about 185 
see . 


Figure 5-17. Tension vs. time for a case similar to Figure 5-16, 
except that the active damping mechanism is turned on at_acquisition. 
Dam£ing parameters: T » 1.6 x 10 dynes. ¥ » 2.66 x 10 cm dynes 

sec p = 0.012 sec*^. 

Figure 5-18. Tension vs. time and radial vs. in-plane behavior 
for the case shown in Figure 5-17, plotted with an expanded time scale 
for comparison with future runs and with an acquisition time of t = 10 
sec. (a) Tension vs. time. (b> A series of radial vs. in-plane plots 
of the tether system. The succ jssive plots are at two-second 
intervals, the same as the plotted points in (a), and the station 
position is arbitrarily shifted to the right at each interval so that 
the horizontal axis is also, effectively, a time axis with the same 
scale as (a). A constant (in time) vector is subtracted from the 
position of each mass to emphasize the displacements. The vectors 
(one for each mass) are snch that at t « 0 the spacing between masses 
is 1/2 km. 

Figure 5-19. Post-acquisition tension variations are reduced by 
turning on the active damper at t » 0 and delaying acquisition until t 
" 9 seconds. 

Figure 5-20. Tension vs. time when the payload is not acquired 
until 20 seconds after the damper is started (at t - 5 seconds). The 
tension had already begun to decrease when acquisition occurred and 
went slack at t « 31 seconds. (Note the expanded time scale in this 
plot . ) 

Figure 5-21. Plots, for the case in Figure 5-18, but with 
acquisition occurring with a 3 meter/sec relative velocity in the 
in-plane direction along the orbit, (a) Tension vs. time, (b) Radial 
vs. in-plane behavior. 
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Figaro S-2?. As for Figure 5-21, except that the velocity is in 
the out-of-orbit plane direction, (a) Tension vs. time, (b) Radial 
vs. in-plane behavior, (c) Radial vs. out-of-plane behavior. 

Figure 5-23. As for Figure 5-21, except that the velocity is 
radial, directed away from the station (downward), (a) Tension vs* 
time, note the tension variation is greater than the other cases, (b) 
Radial vs. in-plane behavior. 

Figure 5-24. As for Figure 5-21, except that the velocity is 
radial, directed toward the station (upward), (a) Tension vs. time. 
Note that the tether segment nearest the payload goes briefly slack. 
This may be an artifact of the lumped-mass model, (b) Radial vs. 
in-plane behavior, the dashed lines indicate slack tether segments. 

Figure 5-25, Co-ordinate relationship in dumbbell configuration 
used in orbital pumping study. 

Figure 5-26, A rough planner sketch of the orientation of 
Mercury's axis of minimum moment of inertia, at different points along 
its orbit, given that the rotational period is two-thirds of the 
orbital period and that this axis is aligned with the Sun-planet 
vector at perihelion. 

Figure 5-27. Orbital transfer from GEO to LEO of a dumbbell OTV. 
Initially, the system may synchronize in a 9:1 resonance due to spin 
orbital coupling as shown. During the decay, or even from the 
beginning if the characteristics of the dumbbell and initial 
conditions are correct, the dumbbell may low order resonances - 8.5:1, 

8:1, 1:1. It will finally synchronize in the 1:1 spin orbital 

coupling in the final LEO orbit where the retrieval to compact 
configuration will be performed. 

Figure 5-28. Vector relationships. 

Figure 5-29. Spacecraft configurations for jrbital pumping. 

Figure 5-30. Eccentricity of the center-of-mass orbit vs, the 
orbital phase (true anomaly/2n) . A critically damped dumbbell 
configuration, with oscillatory frequency tuned to the orbit, is 
allowed to evolve, converting orbital energy to internal (e.g. 
thermal) energy, circularizing the orb^t. The numerical simulation 
was made with SRTHOOR. 

Figure S-31(a). Tether length vs. orbital phase for one 
particular orbit of the study shown in Figure 5-30. (b) The 

derivative of the curve shown in (a), together with the gravity 
gradient forcing function. These are nearly in phase and their 
product (the work done by the orbit on the dumbbell) is nearly 
maximal, (c) A plot similar to (b) for a case with zero time lag in 
tether response. The work is shown as well, and obviously integrates 
to zero over an orbit. 
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1.0 I^n^RooocTION 

The work carried oat under this contract covered a diverse group 
of topics under the general heading of the study of the dynanics of 
spaoebome tethers. SAO has been working in this area for alaost ten 
years under internal and NASA sponsorship. This effort was directed 
at identifying the state-of-the-art in tether raodelling asiong 
participants in the Tethered Satellite Systesi (TSS) Prograa, extending 
the state-of-the-art to nodel the slack tether and study its behavior, 
and studying certain advanced applications of the tether to problems 
in orbital mechanics. 

This report was assembled primarily from the monthly reports 
submitted under this contract, althongh much new material is included, 
particularly in the sections on high-resolution modelling of the 
tether and advanced applications. Section 2.0 summarizes the results 
of this study. 

Section 3.0 encompasses the review of the features and 
applications of the TSS software set. Software status is as reported 
to SAO in March 1983. 

The dynamics of the tether itself has traditionally been studied 
by modelling the tether as a series of Interacting lumped masses 
connected by springs and dampers. The resnlting equations of motion 
are integrated using standard software routines such as the Gear 
integrator used in SKYHOO'C. This approach, althongh mathematically 
rigorous, is computationally slow and of relatively limited 
resolntion, typically 10-20 mass points. Section 4.0 describes the 
successful SAO effort to move beyond these limitations by modelling 
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the slack tether analytically with as many as SO mass points and the 
application of this new model to a study of the behavior of a brr''en 
tether near the Shuttle. 

Lontitudinal oscillations can occur in the tether from a variety 
of causes, inoluding control reel maneuvers and thruster o- orations. 
Proper use of the control reel can damp these oscillations efficiently 
without adding additional hardware to the system. Section 5.1 
describes such a reel control algorithm developed by SAO and gives 
examples of its use, including an example (in Section 5.2) which also 
demonstrates the use of the tether in transferring a heavy paylor.d 
from a low-orbiting Shuttle to a high circular orbit. 

Capture of a low-orbiting payload by a Space Station in high 
circular orbit is described in Sec*^ion 5.3. 

Finally, we show that in some eases it is possible seemingly to 
"pull oneself up by one's own bootstraps." Section 5.4 shows how 
energy transfers within a "du *’ bell"-type spacecraft by cyclical 
reeling operations or gravitational effects on the natural elasticity 
of the connecting tether can circularize the orbit of the spacecraft. 
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2.0 SUMMARY CONCLUSION AND RESULTS 

One-ptge tuamaries of conclusions and results are presented in 
this section for each study topic covered by this report. 

2.1 R eview of Tether Ovnamics Software 

19 different tether models are in use: 7 using angle variables 

12 using a finite element approach. 

Little commonality of hardware and software was found among the 
models although modelling approaches are often similar. 

Angle variable models are used for long term studies, such as 
control law development; finite element models address the dynamics 
of the tether itself. 

Only TE].*HER REEL (MMC) . "Tethered Satellite" (CDC) . and SKYHOOK 
(SAO) model, the dynamic behavior of the entire 
Subsa tel lit^/Tether/Shut t le System, 

Only SKYHOOK models the dynamic behavior cl entire tether system, 
including its electrodynamic and thermal behavior. 

Improvements in software computational efficiency, particularly 
in modelling the slack tether, will be necessary to meet anticipated 
real-time simulation and safety study requirements. 

New software should be developed with transferability in mind, 
although there is no compelling reason all software need be 


transferable • 
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2.2 Hi-Raaolution Study of Tether Behavior Meat Orbiter 

An efficient two-dimensional analytic model of the slack tether 
called **SLACR2'* was developed under this contract with S04> mass point 
resolution. 

SLACR2 includes a realistic deployment boom model and the major 
external forces* the gravity gradient force* air drag* and the 
Coriolis force. It runs about 100 times faster than SKYHOOK under 
similar conditions. 

The tether break study yielded results similar to SKYHOOK runs* 
but the higher (2-3m) resolution provided new insight into ^he 
behavior of a broken tether. 

Deployment boom oscillations add energy and oscillations to 
tether motion after break. 




2.3 Study of Advtaoad Uta> pf the Tether 


2.3.1 Algorithm for Damping Tether Longitudinal Oscillations 

Operation depends only upon measurement of tether tension and 
deployed length. 

''Damper^' fully controls tension oscillations within a few cycles 
after release of a heavy payload using reel-in/reel-out deployment 
maneuver. 

Algorithm is '^robust" with regard to parameter sensitivity. 

Algorithm allows deployment of payload using a pure reel-out 
manenver which limits tether tension to pre-release value. 
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2.3 StndT of Advanced Uaea of the Tether 

2.3.2 Payload Tranafer to Circular Orbit 

A 9,S-ton payload «aa deployed from a 100 ton Shuttle using a 
61-km tether of 3-mm Kevlar and the reel-out maneuver and algorithm 
described in Section 5.1. 

Initial Orbit 

Perigee: 199km 

Apogee: 519 km 

Eccentricity: 0.0237 

Pinal Payload Orbit 

Altitude: 575 km 

Eccentricity: 0.00017 

Pinal Shuttle Orbit 
Perigee: 165 km 

Apogee: 514 km 

Tether tension was always lets than stable deployment tension and 
never aero. 


0 



2.3 Stadv of Advanced P»ei of the Tether 

2.3.3 Payload Acqniaition by a Tether Deployed from a Space 
Station 

A 9.5-ton payload was acquired by a 100-ton Space Station using a 
44-km length of 3-ma Kevlar and the damping algorithm of Section 5.1. 

Initial Space Station Orbit 

Altitude: SOO km Circular 

Initial Payload Orbit 
Perigee: 220 km 

Apogee: 456 km 

Acquisition is stable and relatively insensitive to velocity 
mismatch between payload and tethered teleoperator. 
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2.3 Study fif Advanced Daea of tlie Tether 

2.3,4 Orbital Pnaplng 

Ouabbe 11-type tethered spacecraft placed in elliptical orbits 
rill self-oircnlarize the orbit due to tidal forces acting on the 
spacecraft and or due to deliberate reeling operations on the tether. 

Applications of this pumping technique include 
self-olrcnlarizatlon of an OTV orbit at low altitude and geostationary 
capture of a satellite from a high-eccentricity orbit. 

Energy may actually be recovered as electrical power in 
transferring a spacecraft from LEO to GEO using tides developed in the 
spacecraft system by the Earth’s gravity gradient. 




Pate 17 


3.0 REVIEW OF TETHER DYNAMICS SOFTWARE 

3 .1 Method of Software Review 

SAO solicited data from Tethered Satellite System (TSS) program 
participants on their tether dynamics models, programming langnages, 
computers, displays, and non-standard support software requirements. 

Responses were received from all program participants; Aeritalia 
(AIT), Torino, Italy: Control Dynamics Company (CDC), Hnntsville, AL; 

Johnson Space Center (JSC) .Houston, TX; Martin Marietta Corporation 
(MMC) , Denver, CO; Marrhall Space Flight Center (MSFC) , Huntsville, 
AL; University of Padova, Padova, Italy: and Smithsonian 

Astrophyslcal Observatory (SAO), Cambridge, MA. These were reviewed 
by SAO and the results of the review and summarized in an applications 
matrix and in the text below. We also comment on the issue of 
software compatibility, outline potential tether dynamics study 
topics, and make recommendations for software development. 

3.2 Summary Deacrintions of Tether Dynamics Software 

Appendix I contains the completed questionnaires as received by 
SAO and a complete list of the respondees. Here the features of each 
of the models is summarixed from data provided on the questionnaires. 

3.2.1 Aeritalia Tether Dynamics Software 

Aeritalia runs three tether dynamics programs: ASSETTO, POSI, 

and BALLAST. All model the three-dimensional behavior of the tether; 
ASSETTO adds the three Euler angles of the snbsatellite. All assume 
that the Shuttle has infinite mass and that the tether is straight and 



inelastic with distributed mass. Atmospheric drag is modelled« 

Shuttle out-of-plane delta-V can be applied, and arbitrary control 
laws can be used. Orbits can be circular or low eccentric. Tether 
and subsatellite mass variations can be accommodated. 

POST models the subsatellite and Shuttle as point masses and 
includes a movable boom. The three second-order differential 
equations in the model are numerically integrated by a fourth order 
Runge-Kutta integrator with variable step size or by a 
predictor/corrector integrator. In the case of rate tether control 
laws# the equations related to tether speed and acceleration are 
removed and the remaining equations are integrated. All effects 
related to the elasticity and flexibility of the tether are neglected 
and the Shuttle orbit is unperturbed in the analysis. 

ASSETTO adds to POST the three Euler angles of the subsatellite 
to model its rigid body rotation. It also includes the diagonal 
inertial tensor of the subsatellite and can handle misaligned 
subsatellite thrusters. The subsatellite attitude does not affect the 
position of the tether# an acceptable assumption for subsatellite 
diameters small in relation to the length of the tether. 

BALLAST is a double pendulum model (equivalent to a three-mass 
model) which allows four different physical systems to be simulated: 
double tethers (satellite deployed from the subsatellite), two-section 
tethers, a subsatellite with variable attitude, and a 
tether/subsatel lite system with a moveable boom. Two different 
mathematical models are used. The first allows control of the primary 
tether; the second# control of the secondary tether. The systems are 
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first solved with respect to the second derlvstives then numerically 
integrated. Instantaneous control of only one tether at a time is 
possible. Transverse and longitudinal wave propagation effects were 
neglected. 

3.2.2 Computer Dynamics Company (CDC) Tether Dynamics Software 

COC uses a modal representation of the wire and models the 
Shuttle and snhsatellite as point end masses. The disturbances 
modelled Include: 1) Gravity terms using a non-spherical earth model 

and third-body perturbations; 2) Aerodynamic terms using various 
optional density models. Including the Jacchia model; 3) Solar 
radiation pressure: 4)Tether bending stiffness: and S) Thrust 

maneuvers by either end body. The dyuamics of the tether is modelled 
by p' rtial differential equations which are solved by conversion to 
finite difference equations with a selectable number of nodes. The 
equations are integrated using a 7(8) Runge-Kutta-Pelberg 
variable-step size scheme. Rigid body motion of the Shuttle and/or 
the subsatellite are not currently included. Only open loop thrust 
maneuvers are allowed. No deployment boom or mechanism is modelled at 
present. Modifications to correct these deficiencies are expected to 
be complete by the end of May 1983. 

3.2.3 Johnson Space Center Computer Dynamics Software 

PMG (Plasma Motor-Generator) models the tether as a continuous 
cable with mass and internal moments of inertia. The D1 (PMG-Dl) 
version currently in use is a preliminary version. Two-dimensional 
dynamics are provided using a constant (specified) length tether 
operating under variable I X B forces with constant (specified) values 
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for the magnetic field* gravity* and an arbitrary net force vector 
which is the vector snm of any other force terms such as atmospheric 
drag* radiation pressure* and indnced electric forces. Specified 
terms may be varied from an external program to match an assumed or 
calcnlated reference orbit* to include unique conditions* or to 
include chrust or drag current inputs. The d 3 mamic model is developed 
for application on a hybrid (digital/analog) computer system. The 
current (interim) version is discretized for a finite element solution 
on a digital system using a bacicward Euler integration package by 
Hindmarsh. This version is cumbersome for changes to the input values 
and is restricted to constant length and in-plane (2D) dynamics. 

3.2.4 Martin Marietta Corporation Tether Dynamics Software 

Martin Marietta runs three tether dynamics programs: TSSONE* 

TETHR* and REEL. 

TSSONE models the subsatellite as a point mass connected to the 
Shuttle by a straight* inextensionable tether. The Shuttle and 
deployment boom are considered as point masses constrained to the 
defined orbit plane. Numerical integration is used with prescribed 
initial conditions to develop time histories of motion of the 
subsatellite and other data. The model includes atmospheric drag 
effects* gravity gradient torques* and may include in-line, in-plane, 
and out-of-plane thruster effects. Numerical integration is via the 
Runge technique. The Gauss quadrature approach is used to establish 
atmospheric drag effects. 



TETHR models the three-dimensional motion of the system using a 
lumped mass model ('^bead^ model in KMC terms). The co-ordinate syetev 
rotates with the orbit. Integration is by a relaxation method in 
which mass point positions are propagated independently relative to 
the center-of-gravity of the Shuttle-Tether-Sub satellite system in 
short time steps; they are not simultaneously integrated. This 
approach relies on damping within the tether to achieve stability. If 
the damping is too lorn, the Integration ''blows up," In exchange for 
the damping requirement, relatively fast integration times are 
achieved. Inherent in the relaxation method is that it works better 
with a more elastic tether, worse with a stiff tether. The 
integration time step must be decreased as as stiffness Increases, 

For typical values of stiffness, a fairly large (on the order of 
seconds) time step may be used. Typical run times with six mass 
points on the CDC Cyber 750 computer were 100 times faster than real 
time, TETHR uses a spherical gravity model although a non-spherical 
model could be easily added. It includes an elastic tether 
representation with hysteresis data derived from tether test data. 

Air drag, including angle-of-attaok, can be modelled. It handles 
tether "breaks" but not rotation of the end masses. 

REEL is a development of TETHR. In addition to the features of 
TETHR, this model simulates tether reel-in and reel-out and allows 
input of a control law to regulate the reeling process. This program 
was not fully verified at the time of writing although 
reel-in/reel-out had been successfully achieved. For 
reel-in/reel-out, the mass points are "restrung" at appropriate 
intervals to maintain accuracy. Re-stringing and tether cut are 
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■ataally eiolntive In a given run. 

3.2.S Marihall Space Flight Center Tether Dynamioa Software 

Marehall Space Flight Center presently runs four tethor dynanics 
prograns: TETBBR. TETHBR-PD, SKYHOOK (Foraation Flying), and SKYHOOK 

(Tether Configurations). The two SKYHOOK programs are derivatives of 
the SAO SKYHOOK program discnssed further on in this section of the 
report. 

TETHBR simulates the system with the satellite considered as a 
point mass connected to the Shuttle by a tether which hat matt and 
which is extensible. Disturbances modelled include the spherical 
gravitational field, aerodynamic drag with an exponential rotating 
atmosphere, and longitudinal vibration of the tether. Numerical 
integration of the differential equations is by a fourth-order fixed 
step size Runge-Kutta integrator. The current version does not 
include modelling of a deployment boom nor satellite thrusters. 

Changes to incorporate these features are expected in the near future. 

TETHBR-PD is a two-mass-point model which simulates subsatellite 
trajectories in three dimensions. It assumes a circular orbit and 
incorporates a deployment boom oriented along the local vertical. The 
tether is modelled as a massless inelastic wire: no end mass rotation 

nor atmospheric drag is included in the simulation. Retrieval of the 
snbsatellite is modelled by a control law which includes a tether 
retrieval rate limit, tether angle and angle rate thruster control 
law, and which hat the tension control gain adjusted for short 
tethers. Integration of the equations of motion is by a fixed step 
size trapezoidal integrator. NSFC describes this program as a fast. 
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• impl* tianlator of aatellite t.rajactorlea which it the work horse for 
studying control laws . satellite inpulte reqnireaents, and gross 
systesi characteristics. 

S1CYR00K (Pornation Plying) nses the tame environBental model as 
SAO SKYHOOK which is described below. The tether model was changed to 
permit individnal mast elements (sobsatellites) to be in separate 
orbits with state vectors referenced to a main spacecraft (i.e., a 
Space Station). Additional program changes were made to allow the 
program to run on a UNIVAC cosipnter. 

SKYHOOK (Constellation) is the same as SAO SKYHOOK with minor 
changes to npdate tether system confignration and control laws. It 
too has been adapted to rnn on a UNIVAC machine. 

$. 1.6 University of Padova (Italy) Tether Dynamics Software 

The University of Padova rnns the original version of SAO SKYHOOK 
modified to provide polar coordinate ontpnt and to accept Shnttle 
orbital element variations. It does not inolnde the SAO 
electrodynamics package. 

3.2.7 Smithsonian Astrophysioal Observatory Tether 
Dynamics Software 

SAO runs a variety of general purpose and specialized tether 
d/aamics models. Integrators, data processors, and display programs. 

SKYHOOK models the tether as an N-mass point system where N is at 
most twenty in the present software. A ten mass point model runs in 
approtimately real time. Bach mass point b*s the physical 
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ohartcttrlitlos of ■ finite Ien|th tether segment and mass point 
interaotiona are determined by a set of coupled differential eqnationa 
which are numerically integrated to determine the behavior of the 
system. There are no restrictions on mechanical, thermal, electrical, 
or bulk characteristics of any tether segment nor are there any 
restrictions on the initial conditions representing eccentricity or 
orientation of the Shottle orbit. SKTHOOK models all potential 
physical and dynamical effects on the tether/snbsatellite/Shntth 
system. The namerical integration is done in Cartesian oo-ordina> . 
Rotation of the snbsatellite was added in a recent upgrading of the 
software. 

DUMBEL models the tether system as a massless elastic wire with 
two end masses. It includes position, velocity, and rotation of the 
subsatellite in the integration variables. Numerical integration is 
done in Cartesian co-ordinates. DUMBEL also computes initial state 
vectors as a function of orbital and tether system parameters. Wire 
mass points, and their initial conditions, can also be generated along 
the wire to set up initial conditions for SEYBOOK. 

STABLE models the tether as a rigid massless rod. It assumes 
that the Shuttle is in a circular orbit and integrates the in-plane 
angle, out-of-plane angle, distance to the subsatellite, and the rate 
of change of each as a function of time. The velocity of the 
subsatellite ?sn be replaced by a rate control law. 

WAVES is a "many'* mess point model of the tether which integrates 
only the radial variable with no gravity field and, hence, runs very 
quickly even with large numbers of mass points. It has been run with 



•t mtay at 40 matt poiota and is typically nsad in special case 
stndies to yerlfy or to form a basis for studies with other software. 

SLAC^2 it a mnlti-aass point free~space nodal of the tether for 
studying the behavior of a slack tether with exceptionally high 
resolution and with high conputational speed. Mass points can be 
arbitrarily distributed along the wire and concentrated in areas of 
Interest. As nany as SO mass points have been used in recent studies 
with this node! but expansion to 100 or more nass points appear to be 
straightforward. It runs SO to 100 tines faster than SP'ROoC in 
equivalent slnulations. The nodel assunes that the center of nass of 
the tethered systen is in circular orbit and that the co-ordinate 
systen rotates with it. Gravity gradient, Coriolis, and air drag 
forces are included. It is presently written in two- dinens ions; 
extension to the third dinenslon is not difficult. Forcing equations 
have been linearized (expanded in the first order in the snail 
quantities x/R and y/R where x and y are the co-ordinates of the orbit 
center and R is the orbital radius. 

GRAPH. PLOTFIL. PL0TFIL2, INIT, CRONCH7, and CRUNCH9 are routines 
used with each of the SAO tether dynamics prograns to process and 
display data. 
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3 . 3 Tethor Pyntmiot Software Apolicationa Matrix 

3.3.1 Tether Dynaalcs Software Applications Areas 

Potential applications of the tether dynamics software set are 
listed across the top of the Applications Matrix (Flgnre 3-1) . The 
matrix identifies those application areas considered most suitable by 
SAO for each program described In the questionnaires. Topics of 
Interest within each of these areas are described below. 

Safety Studies: 

The behavior of tether '‘.fter a break, reel Jam, or various 
control system failure modes, as well as the development of methods 
for dealing with such situations, are inclnded in this category. 

Among the results of such studies would be appropriate Shuttle 
avoidance maneuvers, the specification of tether oscillation damping 
algorithms for the control reel, the design of tension wave shock 
dampers for the snbsatellite, and the development of **safe operating 
area" curves for the deployment, retrieval, and station-keeping 
mission phases. Such corves would identify the safe Holts of tether 
velocity and departure angles as a function of orbital altitude, 
tether break strength, tether deployed length, snbsatellite mass, 
snbsatellite thruster activity, STS mission safety requirements, and 
other such applicable parameters and requirements. 
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Otbiter l>rnamic8 ; 

The Shuttle ceater-of-aeee moves off its normal position as the 
tether is deployed, remains so daring station keeping, and returns to 
its normal position as the tether is retrieved. The behavior of the 
Shut tie /Tether system must be aocnrately modelled under eaoh of these 
conditions to determine the effects of Shuttle and subsatellite 
thruster firings and tether control operations on the behavior of the 
combined system and to evaluate contingency maneuvers in tether 
failure situations. 

Sub satellite Dynamics ; 

Detailed study of the motion of the subsatellite while it is 
deployed is essential to baselining its performance as an experiment 
platform. These studies would determine the behavior of the deployed 
subaatellite under dynamic inputs from the Shuttle, tether, 
atmospheric drag, and the action of its own thrusters. Attitude 
control requirements and performance specifications could also be 
determined by means of these studies. The behavior of the 
subsatellite in the deep atmosphere is of particular interest because 
of the high drag level and potentially high thermal load anticipated 
in such a mission. 

Cp.Bltal ItiS Qptimixation ind Verification : 

Control lav optimixatiou will continue to be an area of interest 
as the overall Shnttle/Tether system dynamics are further studied and 
evaluated. This area also includes passive and active retrieval 
studies (retrieval of a subsatellite with a failed in-line thruster is 



one potential study), stabilization of the system during station 
keeping, and deployment and retrieval time minimization. 

Tension wave propsgation and elastic and damping properties of 
tether materials fall in this area of study. Modelling of real tether 
characteristics will be an important part of verifying the true 
dynamics of the tether system. 

3.3.2 Comments on the Dynamic Model Set 

Coaipnter programs for stndying the dynamics of tethered satellite 
systems fall into two main categories. In one, the Shnttle is assumed 
to be in a fixed orbit and the equations of motion of the snbsatellite 
relative to the Shnttle are used to integrate the motion of the 
system. The position of the snbsatellite is specified in a spherical 
coordinate system by two angles and a radius vector. The angles are 
the in-plane and ont-of-plane deviations of the snbsatellite. Since 
the motion of the Shnttle and the dynamics of the tether are not 
considered, these computer programs are compote tionally efficient and 
can be need to stndy the behavior over long time periods. The primary 
application of such programs is the study of deployment, retrieval, 
and station keeping. Control algorithms for the reel motor or 
thrusters can be developed and tested to optimize processes snch as 
the retrieval of the snbsatellite. Rigid body rotations of the 
Shnttle or snbsatellite may also be inelnded in the models to iitndy 
attitude control. This type includes the Aeritalla and MSFC models, 
and TSSONE used at Martin-Marietta. 
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In the other epproech, the tether system is represented by • 
discrete set of masses or segments. The motion of the discrete 
elements is then integrated as a function of time. This approach is 
more general bnt relatively more expensive in terms of computational 
time in proportion to some power (2.0 to 2.S) of the number of 
elements used to represent the system. It has the advantage, however, 
of allowing the orbital motion and dynamics of the tether itself to be 
studied and is the most accurate representation of the tictxinl system. 
Rotation of the Shuttle or subsatellite may also be included to study 
attitude control. A variation of this approach is to use normal modes 
to describe the tether dynamics. The detail achieved depends upon how 
many modes are integrated. Deployment and retrieval can be studied 
with all these models by including control laws and a method of 
handling the variation in tether length. Models of this type are 
SCTHOOK and the CDC tethered satellite program. Program TETHR written 
at ^fMC dispenses with the matrix inversion usually required by the 
integration routines to achieve fast Integration. However, this 
method is stable only if sufficient damping is included in the model, 
although this is only a minor constraint in many modelling sltnations. 

Only TETHR/REEL (MMC), "Tethered Satellite" (CDC), and SRYHOOR 
(SAO) are capable of addressing the dynamic behavior of the tethered 
satellite system as a whole. Only SRYHOOR can integrate the dynamic 
behavior of the entire system with its electrodynamic end thermal 
behavior, including all potential dynamic inputs. 
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3.3.3 Commeata on Hardware and Software Compatibility 

There are few polnta of hardware and aoftware commonality among 
tethered aatelllte dynamiot model nsera and developera (See Figure 
3*2). RBm. (NNC). TETHBR-PD (MSFC) . and SCTHOOt (SAO) were 
speoifically designed for transportability, bnt even these were 
written in two different langnages, ALG0L68 and FORTRAN IV, and run on 
three different machines: CDC Cyber NOS-1, IJNIVAC 1100, and the DEC 

VAX. FORTRAN IV is the most coaimon language used; other versions of 
FORTRAN are next most common. Only JSC and MMC nse the same computer, 
the UNIVAC 1100. There is no commonality among display systems. All 
users can handle magnetic tape for transfer of data and programs, 
however. Even the programs which were designed for transportability 
would require support from the developer to get them running 
elsewhere, since some reprogramming would be required to meet the 
requirements of different computer and display systems. SKTHOOX was 
successfully transferred to MSFC and the University of Padova in this 
manner. A SXTHOOC tape has also been sent to JSC and should be 


running there shortly. 


Figure 3-2. Tether Dynamics Software Computer Facilities and 
Languages in use. 
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4.0 HIGH RESOLUTION STUDT OF TETHER BEHAVIOR NEAR ORBITER 
4.1 Modal Development 

Under » prior contract (NAS8-33591) , atndiee were done on the 
behavior of the tether under various failure modes such as Jamming of 
the reel motor during deployment and a break at some point along the 
tether. Various analytic and numerical techniques were used to study 
the behavior under these extreme conditions. The analytic study 
modelling the tether as a single lump suggested that the recoil 
velocity after a break Is independent of tether length. This result 
was confirmed by numerical Integration under the assumption that the 
recoil is purely elastic with no damping losses. The numerical method 
used the SHYHOOE program which models the three-dimensional orbital 
motion of the tether by breaking it op into a series of discrete 
segments. The discretization procedure makes it possible to model a 
slaok tether, however, the numerical integration becomes much slower. 
In order to obtain good resolution, it is necessary to have as many 
mass points as possible. Unfortunately the amount of computer time 
required rises rapidly as the number of mass points increases. The 
increase is partly due to the higher frequencies that exist when the 
wire is modelled in smaller sections. Damping can be used to 
eliminate the high frequencies and obtain faster integration using 
methods designed for stiff systems. Of course, damping suppresses all 
longitudinal modes along the wire and cannot be used to study 
longitudinal motions if the behavior is nearly elastic. 



Siaot blgh retolatlon and high compotational speed are necessary 
for efficient study of snch failnre modes, alternative methods were 
sought to these standard approaches. For very short pieces of wire, 
the motion after a break is nearly one dimensional and force free, if 
Coriolis and gravity gradient effects are neglected. Therefore, a 
one-dimensional program called WAVES was written under the previous 
contract to allow more economical integration with a larger number of 
mass points. The program waa naed with 40 points and the results were 
presented in the Interim Report to contract NAS8-33591 dated March 
1982. The one-dimenaional program was particularly useful in studying 
stress propagation and recoil velocity effecta as the wire goes slack 
both under purely elastic conditions and with various levels of 
damping. Once the wire goes slack, the points representing the masses 
are in free flight most of the time, coming into tension only when the 
distance between mass points becomes equal to the natural length of 
the wire segments joining the masses. Bach time a segment goes in or 
out of tension it causes a diacontinnity in the numerical integration. 
With a large number of masses the integrator is continually 
encountering discontinuities. The time at which the discontinuity 
occurs must be found by trial and error and the polynomials used I’or 
extrapolation muat b<> rebuilt past the discontinuity. 

The awkward manner in which WAVES dealt with discontinuities led 
us to develop Sk^other approach in which eacL bounce is treated as a 
point event and iti which both the free flight and the bounce are 
handled by analytic techniques. This development led ultimately to 
SLACW2. the model used in this study. 
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Initially, the program was written in one dimension only to 
oompnte the force-free motion of the section of wire next to the 
Shuttle. The program used the same main rontine need in WA\'ES to take 
advantage of the inpnt and ontpat facilities already developed, but 
used a new analytic snbrontine written to have the same main program 
interface as the integrator. The main routine computes the recoil 
velocity analytically from the system parameters and assigns this 
initial velocity to each wire mass point. The subroutine then 
computes the relative velocity between each pair of mass points and 
the time required for the points to separate by a distance equal to 
the natural length of the segment of wire connecting the points. The 
pair of points having the smallest time to the next bounce is 
determined and all the positions updated to that time. The subroutine 
then returns a linear polynomial giving the position of each mass vs. 
time. 


Initial test rnns done with this analytic program showed a large 
percentage of simultaneous bounces, an unexpected and almost certainly 
non-physical result. When two bounces are simultaneous, the order in 
which the bounces are treated in the program can be determined either 
by convention (depending on how the computer code is written) or by 
roundoff error (where numbers cannot be represented exactly by a 
finite number of digits). 

As detailed in Monthly Reports 5 and 6 to this contract, a set of 
rule I was developed for handling these bounces in a way which led to 
results in this model consistent with normal integration techniques. 

A more general two-dimensional formulation of the model was then 
developed which yielded the immediate advantages that the symmetries 



leading to aiamltaneoua bonnoet (and the related dlfficnltlee) in the 
earlier model were broken, that a tether deployment boom offset from 
the local yertical could be included, and that the influence of 
gravity gradient. Coriolis, and air-drag forces was readily modelled 
and interpreted. 

4.2 at SLAC1C2 Mode: 

4.2.1 Model Description 

SLACK2 is a "lumped mass" model, in which the physical tether is 
divided into segments. The properties f>t each segment are mapped onto 
a spring with a mass at the end away from the point of attachment. 

The point of attachment to the Shuttle or tether boom is modelled 
separately. 

The model allows the tether to be divided into unequal length 
segments. It uses a coordinate system rotating with the tether 
system's orbit center and assumes that the orbit center is in a 
circular orbit. It is two-dimensional in the plane of the orbit. 

The prime charaoteristic of this model is that it assumes a 
tether which is slack almost all the time. The springs connecting the 
masses are never extended past their natural length and when a spring 
does become taut the masses on its ends rebound instantaneously. The 
rebound is computed in the center-of-mass system for the two masses as 
a reversal of the momentum coatponents along the line Joining the 
masaes. Dampfng has been included by allowing the velocity to be 
reduced by a constant percentage at each rebound. 
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Tb« teaaon for rettrlotini thl* model to a alack-wlth-bonncea 
eate li computational almpllclty. It leada to eqnatlona of motion for 
the maaaea which are atrlctly decoupled tetween rebounda. The 
aolutlon for the motion of each point then becomea atralghtforward aa 
deacribed In Section 4.2.2, below. 

Initial condltlona In the preaent model are for a tether deployed 
▼ertloally upward, with a uniform recoil velocity toward the Shuttle 
auch that the elaatlc energy of the atretched tether before the break 
la converted to kinetic energy. Simnlationa with the full SKYHOOK 
program deacribed In prevlona reporta ahow that thla la a good 
approximation to the atate at the moment when the broken tether goea 
alack. The program la modular, and theae initial condltlona could 
eaally be generalized. The aegmenta are of unequal length to break 
the aymmetry leading to aimultaneoua bouncea. 

SLACK2 inclndea a flexible deployment boom aa the point of 
attachment. More general boom modela, e.g. damped or rotating, could 
be added eaally. The preaent boom model la purely kinematic, with no 
allowince for the inflnencf of the tether on the boom/Shnttle ayatem. 
The governing eqnatlona for the motion of the boom were derived under 
the aaanmptlon that the boom, which had been bent by the tether and 
waa in equilibrium before the b' ak, la releaaed anddenly at the 
moment Integration beglna and oaelllatea with no damping. Relaxation 
of theae reatrlctlona moat await a more complete tether model allowing 
the computation to proceed from the Initial break. The Shuttle la 
aaaumed to be an Infinite maaa at the origin In the co-ordinate ayatem 
of the model. The boom oharacteriatica were aa provided by MSFC: 
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ltA|th: 849 Inchai, weight: 104.2 lbt.» binding stlffn^sit: (El) 1.3 

z 10® lb. z in^. 

This model alio inolndea the ezternal forces: air-drtg« 

Coriolis, and gravity gradient. These forces, in the approzimations 
used, all result in li. ear differential equations for the coordinates 
of each mass; the masses remain uncoupled. The drag is modelled as a 
balloon attached to each mass point. Thi# choice, rather than the 
SKYHOOK moael of a taut segment between adjacent masses, was made to 
decouple the equations of different masses. For the slack case, it is 
believed to be also physically more realistic. 

For display purposes SLACK2 creates a file in the same format as 
the SKYHOOK plot file so that the ezisting pcst~processor ran be used 
to generate plots in the familiar side-view (radial vs. in-plane) 
format. 

Two distinct though related concerns of utmost importance are the 
robustness of the computational modal and the stability of the actual 
physical aystem (in the sense that minor changes in the initial state 
do not lead to drastically different states at later times). These 
issues are discussed below along with a review of the mathematical 
basis for the model. 

4.2.2 Details of Mass Point Motion 

The forces acting on the mass points which simulate the tether 
are linear approzimations. Coriolis and gravity gradient forces ars 
given by the first-order terms of expansions in z/R and y/R, where z 
and y are coordinates in a Shuttle-centered system and R is the orbit 
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radius. Drag fore* is aodelled using a ooostant ataospherio density 
with a strsaalng velocity auoh greater than the tether's Internal 
velocity. The aooeleratlon produced by drag, as well as by 
gravitational forces, does not depend on the segaect length modelled 
by a given aass la this aodel. The cross-sectional area, hence the 
foroj. Is proportional to the length, and therefore. In proportion to 
the aass. The acceleration dne to drag Is computed as 

3 - (1/2)F C_ Ap v*/a 
u ^a 

where A ■ (diameter) z (length) Is the segment's cross-sectional area: 

F « 0.75 is a factor to account for projection; « ?. is the 

coefficient of drag: Is the atmosphere density; v is the 

streaalng velocity: and a is the mass of the tether segment. We used 

-13 -3 

p^ ” 1.35 z 10 ga CB . the value given by Jacchla (1971) for a 
1000° ezosphere at 220-ka altitude. 

The coordinates z and y are defined relative to the Shuttle (or 
the orbit center of the system if the Shuttle is finite aass). z is 
the outward radial direction, and y Is in the forward along-orbit 
direction. Motion perpendicular to the orbit plane Is decoupled froa 
In-plane motion and is ignored. The equations of motion are then: 

z » 3<l^z +2Qy 

y - -20i-0 (1) 

11/1 

We scale the time by defining r ■ Ot. where <1 “ (GN/k )^'^ is the 
orbital angular velocity and then. 



3* + 27' 


( 2 ) 
(3) 

3 0 0 2 

0 0 0 1 

0-200 

T 

and d “ [0. 0. 0, D] . Equation (3) ia onr fnadaaental equation. 

Tliis ia an in-hoaofeneona linear ayatea of firat-oxder 
differential eqnationa with oonatant coefficients and ia aolnble by 
atandard teohniqnea (e.g., Boyce and DiPriaa. 1977, Chapter 7). In 
outline, one flnda a partionlar aolntton P^(t) to Equation (3), and 
then any other aolntion aay be written 


original 

OF POOR QUALITY 


y - - - 2*'-0 


where ( )' ■ d( )/dt and 0 =■ 0/0^. We aay re-write thia aa 


p' “ 0 p - d 


T 

where p(r) ■ [z. z*, y, y'J , Q ia a conatant aatriz 


0 10 0 


p(t) » Pj(t) + Pjj(t) (4) 

where Pj^(t) solvea the hoaoteneona avatea 

v\ • Q P^ (5) 

Note that aai partionlar aolntion aay be need. By inapectlon. 
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P^(t) - [-2Dt. -20, 1.5 Oc*. 30x1^ (6) 

To find all aolntlont of (5), which Is fonrth-order. we need four 
Independent solutions, p^^ (r) , P2('c), Pj(t), p^(t); then any solution 
■ay be written Pj^(t) ■ a^ p^ (t) + a^ pj (x) + a^ pj (x) + a^ p^ (t), 
with constants a^ depending on the initial conditions. To find the 
p^, attempt solutions of the fora p(r) = r e^^, where the constant X 
and the constant wector r are to be found. This leads to the 
requirement 


a r • X r (7) 

i.e. X and r form an eigenvalue-vector pair of the matrix Q. 

Q has eigenvalues 0, 0, i, and - i: the multiple eigenvalues 

lead to a somewhat more complex form of the corresponding solution 
(see Boyce and DiPrima, Section 7.8) and the complex exponentials 
corresponding to i and - i may be combined to give real solutions 
with sin(T> and cos(t) terms. The end result is: 

P2(‘c) • 10, 0, 1, 01^ 

PjCt) • [0, 0, 1, 01^ t + [-2/3, 0, 0. 11^ 

Pj(t) ■ [0, 1, 2, 01^ oos(t) + [1, 0, 0, -2]^ 

P^(x) • [-1. 0, 0, 2]*^ cos(t) + [0, 1, 2, 01^ (8) 

Combining '4), (6) and (8), we see that the general solution is of the 
form 
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p(t) ■ 
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2 

t ♦ bt + OT + d coi(»") + • tin(t) 


(9) 


where e. b. o. d. e ere veetore deterained by the laitlel eoaditioae 

tad the dreg. There ooeffioieat vectors ere eoapated es follows. 

First, glvea the ialtiel eoaditioae z . a y , y eoapnte e 

0 0 0 0 


vector a (of leagth 5): 


0 

-6 

0 

3 

0 


-2 

0 

1 

0 

0 


I 0 
0 -3 
0 0 
0 2 
C 0 




ead thea use a to get the coefficieats: 


1 
i 0 


0 -2/3 0 

0 0 0 
0 
1 


0 

0 


0 0 
0 -2 


0 

0 


0 

0 


( 10 ) 


(lie) 


i 0 

i 0 

I 

i 0 

i 

0 


0 

0 

1 

0 


0 

0 

0 

0 


0 -2 

0 0 


(11b) 


0 

0 


0 

3 


0 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 


0 

0 


0 

0 


(11c) 


0 3/2 
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aid) 


0 10 0 


0 2 0 0 


0 0 2 0 


0 0 1 0 01 


( 11 «) 


0 0 0 1 0 


e<* |0 002 0 a 


0 0-200 


la practloo it provat ooavanieat to fora ona aatriz witb thata vactora 


as oolnaas: 


U « (a : b : c : d : ^ 


and ooapnta tha stata Taotor as 


p(t) ■ U oos(t) 
sinCt) 

(Nota tbat TJ is aostly fillad, haTiag only foot zaro alaaants.) 


Rontinas which coapnta tha aass point position osing tha abova 


foranlatioa wara tastad by ooaparing a saapla rnn (followad for abont 
20.000 saconds) with a diraot nnaarioal intagration of (2) using a 
fonrth-ordar Rnnga-Rntta aathod. 
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4.3 Mgdtiliag BiuJIi 

4.3.1 Coaparltoa of SLACC2 with SCYHOOf 

For ttodyiag SLACC2 a ease whieh had pravlontly baaa raa atiag 
the SYTHOOC prograa was aaad. Tha reanlta of the SCYHOOC raa are 
ahowa la Figaro 4-1. Roro, 100 ha of 0.2-oa Caviar tathar with a 
300-kg aabaatallita was iaitially daployad apward. It was cat 200 
aatars froa tha Shattla aad tha raa is startad with a aaifora taasioa 
of 12 kg oa tha tathar which qaiokly (0.03 sacoads) loads to a racoil 
valooity of aboat 3.8 a/sac. Tha racoil was followed for 66 sacoads. 

Figaro 4-2 shows a siaalatlon of a siailar casa asiag SLACC2. 
SLA(X2 iaoladas a daployaant booa which is vertical ia this axaapla 
aad haaca does aot oscillate laterally apoa ralaasa. SLACC2 bagias 
tha siaalatioa with tha tether already slack aad racoiliag. Tha 
recoil valooity is osloalated by assasiiag that the stored elastic 
aaergy is oonvartad aaifomly to kiaatlc aaergy, leadlag to 

Vrac ^ (1 + Hj./2M5g)(30^LMgg)/(Av'^) 

where Hj, is tha tathar aass (before tha break). the sabsatallite 

■ass. L tha origiaal length, A the area. p(sl.S) the density, and E 
12 

(*0.7 X 10 ) the elastic aodalns. To reprodnea tha rasalts of Fignre 
4-1. wa were forced to sat Hj, * 0. Ia cosiputlag the initial tension 
for tha SCTHOOC ran, tha siasa of tha tathar was naglactad. 

Fignre 4-3 shows a ran with tha tathar divided into 40 sagaants. 
bat otherwise tha saaa as la Figaro 4-2 (with Hj. * 0) . Tha gross 
behavior is siailar. bat coaparisoa clearly shows that tha 
point-to-point stractara seen in Figaro 4-2 (and prasaaably that ia 
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large numbers of segments* the plotting program was modified to colt 
plotting the **nass points" other than the boon tip and Shuttle. 
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4-’l) Is Istffflf spurious. Tbis situstion liktlp occurs st sll levels, 
end Is probsblf s fssture of tbs bouncing asss aodsl rstbsr tbsn its 
psrtioular rsslisstiou in SLACT2. 

Tbsrs is ovsrsll sgrseasut bstvssn Figures 4-1 end 4-2. tbst is. 
between SCYBOOC end SLACC2. However, there is cousidersble difference 
in detail. The two sianlations should be investigated further to 
deteraine what is causing the difference. Oddly enough, thongh the 
four segaent SLACC2 rnn shonld approxiaate the SnfHOOH run. the 
40-segaeut run exhibits behavior anoh closer to STYHOOC. 

Figure 4-4 shows the results of a run in which the tether aass 
was included in the coaputation of the recoil velocity; the booa is 
deployed 30^ forward of vertical, but otherwise the physical situation 
is the saae as above. Thirty tether segaents were used and the run 
proceeded for 240 seconds. Note the qua si-periodic notion (the tiae 
between peaks is actually decreasing) and the appearance of daaping. 

It is not clear if energy is actually being lost (to air-drag or the 
booa) or if the aotlons are siaply beooaiug less ordered. 

4.3.2 Running Tiae Considerations 

The fonr-segaent run (see Figure 4-2) took 3.3 seconds of (7U 
tiae. Including setup and output. The analogous SKYHOOC run (Figure 
4-1) took several hundred seconds. This tiae is due to the very 
discontinuous nature of the tether forces in when the tether is slack, 
which causes the nuaerlcsl integrator to search for each discontinuity 
and restart once it is found. The advantages of this new epproach to 
running fuudaaentally the saae physical aodel are obvious. The 
40-segaent run (Figure 4-3) took nearly 20 tiaes the physical tiae 
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Figure 4-4. In plane ve. radial behavior of the tether plotted 
every tvo aeconda for a total period of 240 aeoonda after a break for 
the situation aodelled in Figure 4-2 and Figure 4-3 except that tether 
MSS Is Included In oosputlng the Initial recoil velocity and the booo 
Is deployed 30^ forward of the local vertical. The first 52 seconds 
of the run is shown above (a); the entire run is shown on the next 
page (b). Mote the tether oscillations induced by the notion of the 
booo and the overall behavior slnllar to daaped oscillations. 



It) 



ORIGINAL PAGE If 
I ! OF POOR QUALITY 



Figure 4-4 (b) 


pt|* SI 


•odtlltd, b«t taoh high ratolntioa could not be tohitTod by SCYROOC 
uador any conceivable eoaputatloa budget. 

To study the depeudeace of eoaputatloa tine on the auaber of 
aegaenta the aaae phyaleal ease with 20-. 30-. 40- and SO-tether 
aegaeata vaa run. These results are ahoirn la Figure 4-S. The tlalag 
reaulta are plotted la Figure 4-d where (e) shows the CPU tlae divided 
by the physical tlae aodelled and (b). shows thj auaber of bounces per 
physical aecoad and the CPU tlae per boaaee. The ooaputatloaal cost 
laereaees rapidly, about CPU«(# aegaenta)^, bat Figure 4-4(b) 
deaoastrates that alaost all of this lucraese la due to aa Increasing 
auaber of bounces. The tlae to compute each bounce Is a very weak 
fuaetloa of the auaber of segaeats. 

The eoaputatloa tlae depends oa aore than slaply the physical 
tlae aodelled. The dlseretlzatloa used appears to effect the running 
tlae, even for the saae auaber of segaeuts. As the resolution near 
the boon Is aade finer, the oosq|>utatlonal cost uses. 

Another factor which effects the eoaputatloa tine is the initial 
recoil velocity. A series of precisely eor,parable runs Is not 
available but the ratio of CPU to physical tlae seeas roughly 
proportional to this velocity as alght be espected. The tlae taken to 
reach a particular tether eoaf Iguratlon Is roughly Inversely 
proportional to the Initial velocity. The auaber of bounces is a 
function not of a unit of tlae, but of a unit of "action.'* 

It should be noted that the prograa could be aade to run 
aoderately fatter. Efficiency la prograaaer tlae was eaphaslsed la 
developing SLACT2 and there are anaerous places where changes could 
laprove the eosiputatlonal speed. Sone laaedlately obvious ones are 
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Figure 4-5 (e) 
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Figure 4-6. Tlaing results for the four runs shown in Figure 
4-5 « plotted as ousulstive results for the full 38 seconds of each 
run: (a) The ratio of the ooaputation (CPU) tine to the physical time 

■odelled; log-log plot vs. the number of segments; and, (b) The 
ratio of the number of bounoes to the physical time; log-log plot vs. 
the number of segments. Also plotted in (b) is the CPU time taken to 
obmpute the average bounce, which is seen tc increase much more slowly 
with increasing resolution. 







Pag* S9 


not ooapQtiag th« fall ttata vaotor whoa oaly x aad y ara aeadad 
(whioh woald haWa tba eoat of thaaa oparatioaa): optiaixiag tba 

praaaat foolproof aad alapla. bat ralativaly alow, aabroatlaa whiob 
raflaaa root poaitloat: and aalng a aora aopbiatloatad algorltba for 

oaloalatlag tba aarllaat boaaoa wltb booa tip. Catting tba 
ooapntational tiaa by a factor of two wonld allow tba naa of 30% wore 
aata pointa at tba aaaa ooapntational azpanaa. 

4.3.3 Stability of tba Modal Syataa 

Tbia aaotion addraaaaa tba qaaationa of bow wall SLACK2 
rapraaaata tba real tatbar aad bow pradiotably tba real tatbar 
babawaa. 

On tba road to an oparatiag ooapntar prograa we are lad tbrongb a 
aariaa of lavala of abatraetion. Initially, of oonrae. tbara ia tba 
obTttcal awaten Itaalf: tba laagtb of tatbar (with all Ita 

iaq>arfaotiona introdacad in aanafactara or handling), tba oomplez 
anb-ayataaa attaobad to aitbar and. and tba fall ooapleaant of 
eztarnal foroaa. Aaaaaptlona and approxiaatlona (aneb aa uniforaity. 
tba lawa of aolid aaebanioa. and linaar foroa aodala) ara aada and a 
aatbaaatlcal aodal ia written. For ayataaa anob aa tba tatbar. tbia 
aodal oonaiata of ooaplad aata of differential aqnationa (DE's): 
partial diffarantial aqnationa (PDB*a) for oontinaooa componenta. and 
ordinary diffarantial aqaationa (ODE* a) for diacrata alenaata. Tbaaa 
aqnationa tjrpioally eannot be aolvad exactly and one aaat approxlaata 
tba eolation. One conaon aetbod. partionlarly for eoaplax ayataaa, is 
to ganarata a anatlallT diaoratixad aodal wbleb will oonaiat oaly of 
ODE* a. Tbia aay be dona, for inatanoa, by oonaidaring the ataia of 



Pag* 60 


the oontlsnoaa eomponeata only at a set of dlaorete polnta aad 
approzlaatiag the apatlal derivatives as differeaoes betweea valaes at 
these poiata. The prooedares for (approziswte) solatioas of ODE' a are 
veil eaoagh aaderstood that oa* aov geaerally proceeds to the 
aa«erical ■odel aad its realizatioa as a coswater oroirasi . 

The laaroed mass siodel ased here differs somewhat from this 
approach ia that the spatial discretizatioa is realized directly oa a 
physical level aad we aeed aever write dowa the PDE's of ooatiauam 
meehaaioB. The mathematical model of the lumped mass physical system 
ooasists eatirely of ODE's which are solved as before. 

These stages are aot. of coarse, aeatly segmeuted bat iateract 
oae with aaother. Posseasioa of a well-aaalyzed mathematical model 
will allow oae to determiae more clearly which forces aad effects may 
be aeglected ia prooeediag from the physical system. Compotatioaal 
difficulties ia programmiag the lumped mass model may suggest better 
ways to effect the origiaal physical approzimatioa. 

Whea apeakiag of the stability of the system it is aecessary to 
be specific about what is meaat. The traaaitioa from oae modelliag 
level to aaother eaa either iatroduoe artificial stability (as when a 
critical force or degree of freedom is aot modeled) or coavert a 
stable system into aa uastable oae (damping might be improperly 
neglected, or improper discretization can lead to numerical 
iastabllitles) . Thus, the results of computer simulations must be 
reviewed carefully before conclusions can be drawn about a physical 
tether ia space. 
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la th« laapad ■••• aodal thare are two ebiof araas of eoaoara: 

(1) How pradlotabla la tha bahavlor of tba actual, phyaloal ayataa? 
That la. If we were to conduct actual experiments, would the state 
of tba tathar aftar soma raasoaabla tlaa (say aftar tha first rabonad) 
dapaad atroagiy on tha Initial oonditloaa or variations in tha tathar 
propartiaa? (2) How wall doaa tha (physioal) Inapad swss aodal (point 
aasaaa tiad togathar by aasslaaa tathar sagaants) rapraaant tha 
original ayataa? 

latoitivaly. it would not ba nnraaaoaabla to find that tha ayataa 
itsalf aay show soaa instability. Exparianoa shows that tha 
oonsaqnanoas of a suddan ralaasa of strasa are nnpradiotabla in detail 
and that tha noraal rasponaa of a cut string is to oruapla and move 
about. Howavar. it is not clear whether this represents an 
instability in tha d 3 maaics of tha string or strong variability in tha 
details of tha cutting process. 

With regard to question (2) above, multiple computer runs ware 
doaa to explore tha stability of tha results with raapaot to 
variations in tha details of tha modal used. If somewhat different 
models give closely similar raanlts. than coafidanoa in tha modal is 
high. If models varying only slightly in detail give wildly varying 
rasnlts. tha modal Is in soaa sense suspect. SLACC2 was written so 
that aftar generating a particular discratixatioa of tha tathar, tha 
individual sagmaat lengths could ba varied by (1-t-s) where s is a 
random number in .''ba range, [-0.02. 'fO.021 (2% variation); the 
lengths are than normalised to retain the original total tathar 
length. 
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B«fot 0 •laaittiag tb«M rtaaltt* r«f«r to Flgor* 4-5. which chowt 
tlaulatioac for the teat phyeieel etee aodeled with different noaberc 
(10 to SO) of tether eegaente. Looking nt the final oonflgnretlon for 
etch run (t«3d eeo). we note that though the 10- and 20-aegaent aodela 
bear little reseablanoe to the othera, the 30-, 40-. and 50-eegaent 
aodala yield alallar reanlta. Indeed appear to be converging, an 
Indication of aodel etablllty. 

A 30-aegaent nodal was therefore chosen as a snltable cate for 
the randonlzed partition atndy. Thirty segaents gave about the aaae 
results as 40 or SO la Figure 4-S and ran substantially faster. These 
rnas were oarried through to t«60 seconds to generate sons large-scale 
structure (tether trailing and level with the boon tip). Figure 4-7 
shows the ooaf ignratlon at t«40 and t>d0 seconds for three runs with 
2% segnent length randoa last ion. It should be enphaslaed that these 
are Intended to represent the sane physical tether. 

Figure 4-7(a) shows the configuration at a aonent when the tether 
has extended to nearly Its full length below the Shnttle and rebounded 
for a few seconds: Figure 4-7(b) shows the tether after it has 

rebounded and extended again in its original direction (see Figure 
4-8(a)I. In (a), the consistency between the three results Is 
striking: note particularly that the tips In all three slnnlations 

are beginning to swing upwards in the sane direction. The situation 
In (b) does not seea so consistent, but we anst bear three things in 
aind: the tether extension Is auch less than In (a), l.e. the tether 

Is aore cruapled: the scale Is larger, eaphaslxing the differences: 

and. as a aatter of general experience, all aodels and nuaerical 
processes subject to Initial perturbations will diverge farther froa 
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wave. This saae distinction is seen between Figures 4-5 (a) and (b) 
and Figures 4-5 (c), (d), and (e). Tbe finer resolution results ar 
■ore consistent aong tbe&sel»es and probably closer to the real 
behavior of the tether. 
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tha rafaraaoa aodal aa tia* iaoraaaat ualata tbar* at* tpaeial 
oireaaataaeat. 

la suBaary. if wa do act atpaot datailed aaoaracy of tha aodal oa 
a polat**b]r>poiat batia (aoaathlag which aaaaa ualikaly eoaaidariag tha 
diaerata boaaca aatnra of tha aodal). tha rapaatability of tha 
aiaalatioaa ia aaeoeragiag, 

Thaaa plot^ ooataia aora apparaat datail thaa ia atriotly 
Jaatifiad. Ia Pigara 4-7 (b). thara ara 200 aatera of tathar ernaplad 
iato a apaoa baraly 70 aatata aoroaa. Tha raaolatloa (aagaaat laagth) 
rariaa froa S to 9.S aatara aad ia fiaar aaar tha attaehaaat point. 

Tha diaerapaaoia* ara tharafora oaly a faw aagaaat laagtha at aoat. 
Fiaally, tha atraight aagaaata plottad ara aialaadiag. Tha length of 
tathar batwaaa two aaaa poiata aay ba aavaral tiaaa thair plottad 
aaparatioa aad will ba ernaplad soaawhara in aa allipaa with tha aaaa 
poiata aa fooi. A aora praaiaa raaditioa of tha phyaioal tathar woald 
plot tha ropa batwaaa aaaa poiata ia aoaa way aneh aa a "haxy elond" 
rapraacatiag the naoartaia looatioa of the laagth of slack tathai 
batwaaa tha aaaa poiata. /nat aa ia SKYBOOf. SLAC1C2 ahows where oaa 
aad of each tathar aagaaat is, not the actual tether shape between 
aaaa poiata. Rare, of eoaraa. tha aetaal raaolotioa ia aneh higher. 
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4.4 Conolualoat tad Reoommeadatlont 

If the physioel tystam itself is aastsble thea this form of 
probsbilistio cloud meatioaed above may be the best predictioa oae caa 
make about the tether state aome time after a break. The results of 
these oomparisoa studies are varied eaough that such a coaclusion is 
suggested. Probabilistic results, chough not as clear cut as 
determialstic solutioas. will be about as useful In chat, a range of 
poteatial behaviors caa be derived from the model. This will set 
limits oa the behavior of the real tether; after all, the tether Is 
iahereatly probabilistic, siace it is iahomogeaous and subject to 
non-determiaistic forces such as air drag, random break forces, and 
random attitude control forces. Further study will be necessary to 
determine just what model conditions automatically yield the most 
probable behavior of the tether, and what model conditions are 
appropriate to identifying, say 3a or 5a behivior variations. 

Another corclusion of interest is that this model, at least in 
its present form, is sensitive to input parameter variations. The 
occasional pair of near-aimultaneous bounces is usually handled 
adequately, although there have been cases in which the program gets 
csught in a tight cycle of bonaces (which seem genuine, given the 
model, aid not an artifact of the program). If the entire tether 
comes into tension at the same time, which may easily happen because 
of the linearized forces, then multiple simultaneous bounces are 
possible and cannot be handled adequately in the present model. 

Further thought most be given to this behavior to ensure it is 
correctly handled in future simulations. 



A aajor ftatntt of tha aodal bahavlor, and on# that likely 
appllea to all leaped aa*a aodelt, it the ezlsteooe of two dlatiaet 
aodele. I& high-reaolotioa aodela (large anabera of tether eegaeata). 
•ooa after the elaalatloa atarte a wave propagatea oatward froa the 
aear-Shnttle regloa of the tether reaaltlag la laereaaed velocity of 
the baae ead. la lower reaolatloa aodela. thia effect doea aot appear 
aad the overall behavior of the tether ia decidedly differeat. 
reaaltiag la a acre craavled. "chaotic” coaf Igaratloa. 

Coatiderable work reaaiaa to be doae before thia wave pheaoaeaoa 
ia felly deaerlbed aad clearly eaderatood. Bzperiaeata with oer 
prograa ahoeld deliaeate the depeadeace of the eotoff betweea the 
aodela (aad aech propertiea aa the wave propagatioa velocity aad the 
eaergy traaafer) epoa faotora aech aa the aeaber of eegaeata, the 
diacretixatloa detaila, the recoil velocity, the tether leagth, aad 
other tether proaertlea, each aa deaaity aad diaaeter. Other 
qeeatloaa to be coaaidered are: Why doea the velocity eschaaged 

beooae directed ilfiai the tether? Are there other wavea beyoad thia 
iaitial oae, perhtpa aore aabtle ia their effect? 

Ia aeaaary: 

~The tree reaoletioa of SLACK2 appeara to be a few aegaeat 
leagtha. 

"The tether eoaf igeratioa after oae recoil depeada aoaehow oa the 
way ia which the tether ia dlacretized. However, for large aeabera of 
aaaaea, the flaal coaf igeratioa aeeaa to coaverge. 
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aodel. tad probably tay laaptd attt model, fails to fiad the 
pbysleally sigaifioaat wave pheaoaeaoa for fewer thaa 30 to 50 
aegaeatt. Models with fewer segaeats also show sabstaatlally 
differeat gross tether aotioa aad mask this pheaoaeaa. 

-It is diffioalt at this poiat ia developaeat of the model to 
draw coaolasioas aboat the stability of the physical systea bat it 
appears that at least a probabilistic descriptioa of the behavior of 
the slack tether is attaiaable. 

-The program is ooapotatioaally efficieat. raaniag aboat 100 
tiaes faster thaa SCTHOOt aader similar coaditioas. 

SLACC2 presently models aotioa withia the orbital plane only. It 
woald be simple to add aotioa ia the third dimeaslon aad shoald not 
increase cosipatatioaal time sabstantially. This woald be asefal to 
study (a) instability to oat-of-plane pertarbatioas and (b) effects of 
a sideways-deployed boom (which might be effective la avoidlag tether 
recoil onto the Shattle). 

Farther development woald also allow the springs to cone into 
tension aad remain taat. Since the extensibility of the tether is 
ssMll. the spring forces may also be liaearised. The model woald 
still have a systea of linesr differential eqaatlons (with constant 
coefficients) which may be solved exactly under this condition. 
However, in contrast to the present model, the messes sttsched to taut 
springs woald generate coupled equations. If s segment is else? it 
effectively decouples all the masses above and below it. and these 
problems may be solved separately: but. there is still the potential 

for 4N coupled eqaatioas. where N is the number of tether segments 
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(■«••••). Th« tolatioa of the ooopled eyetea it equivalent to solving 
an eigenvalne/eigenveotor prohlea for the ooeffioient aatriz. Due to 
the fora of the ooeffioient aatriz. SAO believes that a routine oaa he 
developed to eoapute these ezaot solutions unoh aore effioieutly than 
eurrent library routines allow. 

Self-interseotion of the tether is ignored in this aodel. 
Currently, if two sections of tethers approach they siaply pass 
through each other without interacting, a reqnireaent of the 
two-dinensio&«' aodel. Methods of dealing with self- intersect ion are 
under consideration for inclusion in a future three-diaensional aodel. 

Oevelopaent of SLACK2 along these lines is ezpected to lead to a 
ooaputationally efficient high resolution (~S0 aass point) and 
coaplete dynaaic aodel of the tether suitable for a wide variety of 
sianlations and studies. Developaent of the aodel froa this point 
wonld include modelling the dirnaaies of the booa in reaction to the 
tether, the influence of a non-circular reference orbit, a finite aass 
shuttle, and a drag aodel in which the balloon rsJins attached to each 
mass is re-catonlated after each bounce with a projected r-'ornent area 
aodel (this would still leave the equations deconpled betwv 


bounces) . 
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5.0 STUDY OF ADVANCED USES OP THE TETHER 

5.1 A^gRTltla iRs BimUi LgaiilRAiRAl PistlUtigRi ia l&R XcJlRi 

An nlgoritha to control the tetbor reel on the Shuttle bee been 
deelted to dnap out longltndinnl oeeilletlone of the eyeteu. Tbie 
nlgoritba depends only on the force (tension) perceived by the reel 
■ecbnnisa end on the deployed length; it specifies changes to be 
applied to the deployed length. Paraaeters to the control lav reqnire 
knovledge of the tether properties and snbsatellite aass. but no 
detailed knowledge of the systea state is required. 

The fundaaental seheae is siaply: 

- When the perceived tension is high relative 

to soae noaiual tension, reel out aore tether. 

- When the tension is low. reel in. 

This soheae acts to relieve stress on the tether under high 
tension and to add stress under low tension. It will also reaove 
energy froa any tether oscillations, since a full cycle bringing the 
length back to soae original value will cause the tether to do net 
work on the reel aotor. 

Natbeaatioally. we define an applied control length A(t) such 
that the natural length at tiae t is given by 

a^(t) • ♦ A(t) 

0 o 

where is soae noainal or initial natural length and A(t) is given 

A(t) - ftT(t) - T^I - A A(t). 

0 


^7 



T(t) « P«ro«lv«d ttatloa oa t1i« reel at tiae t. 

■ a aoaiaal or tatfat taaaloa: this is a paraaatar 
ohosaa. Idaally. to ba tba aqallibrlaa taasioa. 

C. 3 • ooatrol paraaatars. 

Tha first tara represaats tba fandaaaatal ooatrol seheae of realiag ia 
or oat ia raspoasa to taasioa obaapas. Tba saooad tara is raqairad to 
aaka tba algoritba stabla. Ia tba absanoa of taasioa variatioas it 
aoald oaasa tha ooatrol laagtb to daoay azpoaeatially to aero, while 
if tba aqailibrioa taasioa is iaoorreotly astiaated. it avoids a 
saoalar tara ia tba tatbar laagtb. vbiob woald otharvisa iaoraasa or 
daoraasa witboat boaad. 

Althoagb ia tba oases studied, tba paraaatars T, 3. sad T ware 

o 

bald coastaat (with tha ooatrol law baleg turaad oa sharply), there is 
ao & ariori raasoa why this aust ba so. For iastaaoa. tha ooatrol law 
aigbt ba phased ia aora gradually to avoid a sharp dlsturbaaoa or oaa 
aight vary tba paraaatars as tha systaa properties vary, as dariag 
daployaaat or retrieval. 

A detailed (thaoratioal) aaalysis of this ooatrol algoritba was 
aada whaa it operates on a siapla systaa oonsisting of a aass (tba 
subsatallita) ooaaaotad to a fixed poiat (an "infinite aass" shuttle) 
by a aasslass spring in a oonstant gravity field. This study provided 
a valuable insight into tha systaa behavior as a funotioa of tha 
paraaatars T and 3* A nuaber of factors are left out of this 
analysis: tether aass, partioularly tha finite propagation velocities 

this iaduces; tha fact that tha gravity gradient force does not yield 
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a ooaataat gravity field: Coriolis forces: the Shuttle's fialte 

■ass; aad air drag, eleotrodyaaaio forces, aad light pressure. Of 
these oaissioas, oaly the first seems likely to have aay aigaificaat 
qnalitative effect oa the system behavior, at least ia the limit of 
small oseillatioas. 

5.1.1 Results of the Idealized Case Study 

The model coasists of a spring of natural length z^ and spring 
constant k, fixed at one end with a mass ■ at the other. It is in a 
nnlform gravity field, g. The damping law given above is applied and 
we assume Book's law gives the tension T(t). Let u « /TkTmT, and use 
a displacement L(t) ■ z(t) > (z^ + g/u ), as the dependent variable. 
The equilibrium tension is simply gm, but let ns introduce some 

error s, so that we use T^ ■> g m -t- k s. We then have a system 

** i 2 

L -f u L •• <it A 

A - C ktL- A-sl -BA 

After normalisation of the tine by m we obtain a linear system 
whose characteristic equation is 

k* + (f + B) X* + X + B - 0 

with the dimensionless parameters: 

C • K k/u 

B ■ B/« 

The time has been scaled by w, so the frequencies of the damped system 
are «iX. Since the solutions of the homogeneous part of the system are 
sums of terns ezp(wXt), stability requires that all roots of the 
characteristic equation have negative real parts. So long as B 0, 
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as lahomoieaeous ooastaat t*ra (due to aa exteraal force or to aa 

erroaeoaa estimate of the eqnilibriam teasioa T ) leads oaly to a 

o 

eoastaat term la the solatioa. The behavior of the roots as fnaetloas 
of C aad 3 was ezaaiaed by solvlag the characteristic eqnatioa 
aomerlcally for aaay oomblaatioas of the parameters aad is qaite 
complex (see Figure S-1). The precise quaatity of iaterest- the least 
aegative real part of the set of roots- behaves more predictably: If 

we fix K aad coasider this stability determlaiag real part (let as 
call it simply "the stability") as a faactioa of 3i the stability 
lacreases from 0 at 3 * 0 aad has a cusp at some optimal value, after 
which it slowly decreases. The easp is much sharper for 3 below the 
optimal value thaa above. We aot oaly wish to obtala rapid damping, 
which would lead us to choose the optimal 3, but we wish a "robust" 
damping mechanism, one insensitive to errors in the parameters or in 
the model used for analysis; thus, it is safer to choose 3 to give a 

MS 

damping 3 somewhat larger than optimal, since if the optimal value was 
aimed for but the actual 3 achieved were lower than optimal the 
stability would be seriously effected. More detailed examination of 
the roots suggests reasonable parameter values of C » 1.3, 3 = 0.24; 
for true X in the range 1.0 to 1.5, the damping will not be seriously 
degraded from optimal. 

5.1.2 Numerical Simulation 

The reel control law was implemented in the SKYHOOK program by 
concatenating the differential equation for the deployed (control) 
length with the differential equations for the positions and 
velocities. The control law parameters are read from the input file 
describing initial conditions, and the control length was then solved 
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Flgur« 5-1* Behavior of the roots of the eharact eristic equation 
for fixed values of K. Mote that we plot the negative of the real 
part of each root (pairs of complex roots are Indicated by a dashed 
line) and that for K* $ > 0 all roots have negative real parts (i.e, 
the system is stable).. The quantity of greatest Interest Is the lower 
line at eaoh value of B* This root governs the solution which decays 
least rapidly, and we want to choose paroneters K, B to optimize this 
worst solution. In choosing actual parameters K, B ve may not achieve 
the desired dimensionless parameters, so there Is the additional 
constraint of robustness; slight errors In the paruseters should not 
lead to large decreases In performanoe. 
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•loBg with th« other e^natlOBs otlBg the stiffly stable Oear 
iBtegrator. The goal was to choose dampiag parasieters which led to s 
robBst daapiBg seheae (ooe that behaves well 1b asBy eiroaastaBees if 
BOt optiaally for the Ideally expected case). 

Figure 5-2 shows tSBSiOB variatiOBS whsB a siBusoidal reel 
aaBOBver of period lOS sec is applied with release at 85 see of a 
10-tOB payload, iBitially reeliBg out the tether. This has the 
advaBtage over aB iBitial reel-iB aaaenver in. that the teasioB never 
iaereases over the pre-release equilibrina value: however, experiaent 

with the paraaeters of the sisqile sinusoidal aaneuver being used (see 
the FiBal Report, March 1983, Contract NAS8-33691) failed to prevent 
the tether going slack. Figure 5-3 shows the saae ease where the 
active daaping control law was initiated at t • 100 sec, at which tiae 
the sinusoidal aaneuver was teraiaated. This avoided the slack 
tether, although tension did get close to sero at one point, and can 
be aeea to be daaping the tether oscillations. Froa the sharp turn in 
the tether tension curve when the daaping started, one suspects that a 
gradual initiation aight be beneficial and later runs bore this 
observation out. 

All siaulations used an 80 ka initial deployed length of 0.3-ca 

% —12 2 

Revlar (density 1.5 g/ca . B « 0.7 x 10 dyne/ca ). The aass of the 
shuttle was taken as 100 aetric tons, that of the payload as 9.5 tons, 
and the reaaining teleoperator as 0.5 tons: the tether aass it 0.85 

tons. The systea was in a circular orbit of height 200 ka. 
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Flgurs 5”2. Tension vs, tins for a release study (100~ton 
Shuttle, 9.5-ton payload, 0,5-ton teleoperator, and an 80-ka length of 
0*3-ob Kevlar tether) in which the tether is reeled out first to 
relieve tension. The maneuver is given by (deployed length) - 80 km + 
AIl-eos(2jit/P)l, with amplitude A - -«50 meters, period P - 209 sec. 
After 1/2 period, the deployed length is held oonstant. The payload 
is released at 85 seconds. 
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Figure 5-3. Tension vs. time for a release study In which all 
parameters are the same as In Figure 5-2 except that at t 100 seo 
the sinusoidal maneuver Is terminated (held constant) and the active 
damping control law Is Initiated (parameters T - l.g x lO” dynes; K 
■ 1,6 X lO"* cm-seo"^ dynes“*; - 0.041 sec”*^7. 


I 

{ 

\ 




Pag* 10 


Th« flzad prt-rtlaaaa r««l aanauvar spaeifiad aa addltioa to tha 
aatural langth of Aala(2n t/P i) of atiplitada A and parlod P; tha 
phaaa d was. la tha eataa atndlad. - n/2 ao that tha aaaaaTar atartad 
gradoanp. A tiaa la also apaclfiad at which tha aaBaovar tatalaataa; 
aad, aftar this tiaa, tha aatnral laagth raaalat ooastaat (azeapt aa 
audifiad by tha daaplag eoatrol balow). Oaaaraily, thio atop tiaa ia 
ohoaan aa P/2, ao that tha aaaaavar *'taila off gradually aa wall. 

Tha aotiva daaplag control law ia diffaraat in concept froa tha 
pra-ralaaaa aaaaowar: wharaaa tha aaaauvar ia atrlotly data&aiaad by 

tha paraaatara A, P,d and tha atop tiaa, tha active daaplag algorltha 
raacta to tha atata of tha tathar Itaalf, aa wall aa to paraaatara 
iapat at tha atart of tha aiaulatloa: tha real control laagth 

spacifiad by tha daatpar of aacaaaity dapanda oa tha tathar atata alaca 
tha purpoaa ia to aodify whatever atata ia fonad ia tha diractioa of 
aqallibrlua. 

Wa firat applied tha daaplag achaaa to aaaauvara ia which tha 
tathar ia reeled in before ralaaaa (with a poaaibla aubaaqaaat 
raal-ont), thaa wa iavaatlgatad tha poaaibility of alaply raaliag the 
tathar oat to relieve taaaioa. Tha taaaioa oa tha tathar (at 
aqailibrina) will, of coaraa, be higher with tha 9.5 ton payload atlll 
attached thaa with oaly tha O.S*toa telaoparator (and 0.85 ton of 
tathar). To avoid a raboaad of tha tathar it it thus iaparativa to 
reduce tha taaaioa to aaar tha fiaal daairad aqallibrlua valaa before 
ralaaalag tha payload. Tha raal-ia aaaaavar haa tha advaataga that by 
aaaipalatiag the paraaatara oaa can alto have tha ayataa nearly at 
raat at tha aonaat of ralaaaa. Thia laavat tha poat<>ralaaaa ayataa 
cloaa to aqallibrlua. Oa tha other head, by iait tally pall lag oa tha 
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tcthtr while the heavy peytoed is still sttsched, the nsxiaua tension 
experienced during the nsneuver is Increased over the eqalllbrlnsi by 
about S0%. requiring a stronger tether. If we reel oat rather than 
la. the tension never Increases above the eqnlllbrlua valne; however, 
experlaent with the slaqple deteralnlstle nsneuver always found the 
tether going slack, seeniagly because the systen could not be left 
with a asar-xero velocity and oaclllated strongly after release. 

Our studies of the active danper In oonblaatioa with a reel-ln 
pre-release nsneuver used as a starting point the case shown in Figure 
4-5 of the NAS8-33691 Final Report (Colonbo. 1982). The release and 
sinusoid paraneters for this case had not been optinised. and the 
post-release tension variations are shown in Figure S-4 on an expanded 
scale. (Note that this requires restricting the graph to tines 
greater than 105 seconds after (.he nsneuver starts.) No attenpt was 
nade to Inprove the release t .ne or nsneuver. The active danping 
algorlthn was initiated at the noneat of release, and several choices 
of the danping paraneters are shown In Figure S-S. Figure 5-6, and 
Figure 5-7. Figure 5-5 shows the best danping achieved. Although no 
effort was nade to sea«‘ch fully pataneter space for an optinal 
solution, the behavior Is excellent, with a nlninun tension of about 
16% of the equillbriun value and awjor variations dancing out within 
20 or 30 seconds (conpare to the natural oscillation period of about 
22 seconds in Tigure S-4). Figures 5-4 and 5-7 illustrate how robust 
the algorlthn is: even if the paraneters are chosen very poorly, the 

danping is still effective. In Figure 5-4 the paraneters K and 8 are 
1.75 greater than la the best case Figure 5-5 (this was an attenpt to 
use the neasured oscillation frequency of the systen rather than a 



y 


o 



Figure 5-4. Tension oscillations after payload release, with no 
damping, A ••oel-ln/reel-out sine-law maneuver was used prior to 
re. 3 Bse with amplitude 176 meters, period 114 seconds, phase offset 
-n/2. The maneuver terminated at 114 seconds and release was ;t 117 
seconds. We plot only from t 105 seconds. System: 80~km tether of 

0,3-om Kevlar; 100-metrlc-ton Shuttle; 9.5-ton payload; 0. 5-ton 
remaining teleoperator. 
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spriag/Bast dtrivad frtqaeaey la ooBpatlag daslrad paruietert) . 
Th« dlff«reae« la th« taatloa plots is bsroly disooraablo. thoagh 
Figaro S~S is aarglaally bottor. la Figaro 5-d tho aoaiaol toasioa 
woo oaly <3% of tbo ootaol oqailibrlaa toasioa at tbo Sbattlo; thoagb 
tboro is aa laitial ozoaroioa toward low toasioa, tho bohawior is 
still qaito aocoptablo. 

For tho stadlos of a rool-oat pro-roloaso aaaoavor, wo 
arbitrarily ohoso a poriod P for tho slaasoid of P • 200 soo. (sad 
stoppod tho aaaeavor at t • 100 soo) . Varylag this paraaotor woaid 
probably allow iaproToaoat of tho post-roloase bohavior, bat did not 
proTO aoeessary. Olwea P, siaalstions woro aado for varlons choioos 
of tho aaplitado A aad tho toasioa plottod. Notiag that tho toasioa 
ohaago is olosoly llaoar la A, wo oxtrapolatod to aa aaplitado (A ■ 
>40d aotors) saeh that tho aiaiaaa toasioa (at tho sabsatollito) 
roaohod was tho dosirod post^roloaso toasioa. Tho payload was 
roloasod at tho tlao of this alaiana toasioa (t « 70 soe) aad tho 
daapiag algoritha laitlatod. Tho siaasoid aaaoavor coatiaoed aatil 
t ■ 100 soe to avoid an abrnpt chaego in rool voloeity. Daapor 
paraaotorc wore tho saao as for tho optiaal oaso (Figaro 5-5) 
doaeribod abovo. Tho rosaltiag toasioa is shown in Figaro 5-8 which 
plots tho ooaploto aaaoavor. Althoagh the post-roloaso variations aro 
not so laprossivoly daapod as in Figaro 5-5, the behavior is 
aoooptablo aad doabtloss ooald bo iaprovod with a bottor ehoieo of 
paraaotors (particalarly P) . This oloarly deaoastrates tho 
foasibility of a rool-oat aanonvor eoahinod with an aotivo daapiag 
sohoao. Note: 1) Tho toasioa aovor goes abovo tho pro-roleaso 

oqailibrlaa of 300 x 10^ dyats; this shoald allow a thianor. lighter 
tothor. 2) Tho toasioa does not fall sigai/ioaatly below tho 
oqallibriaa valno. 
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3) Tht tottl tla* for th* Boaeuyor. to a aoar aqallibrlna state. Is 
about ISO see. the saae as for the reel-ia aaaeoyer. 

For the cases plotted ia Figures S-S aad S>8. the best cases, we 
eoaputed the deployed leugth aad reel velocity due to both the fixed 
(slausoid) reel aaoeuver aad the active daaplag ooatrol. For the case 
showa la Figure S-b. theae are plotted la Figure S-*9 aad Figure 5-10. 
The aaiiaua reel velocity la largely deteraiaed by the fixed reel 
aaaenver. the daaplag ooatrol velocity belag relatively saaller aad 
oaly takiag effect vhea the fixed aaaeuver has slowed dowa. The 
aaxlana velocity is 12.8 aeter/seo. For the case showa la Figure 5-5. 
these results were aot plotted aor eoaputed la detail, but the reel 
aaaeuver has aaxiaua velocity 4.9 aeter/see. while the daaplag coatrol 
veloi ity is saall. << 1 aeter/seo. These velocities seea reasoaable. 
but oae could probably decrease thea by varyiag the period of the reel 
aaaeuver. ladeed. If the teasioa reduotioa is roughly proportloaal to 
the aaxiaua aooeleratioa of the aaaenver. we can reasoa as follows: 
the aaxiaua velocity of a aaaeuver A sla(2ff t/P) Is V ■ 2nA/P. and the 
aaxiaua acceleration Is 
a - 4 A/P*, 

or 

V - P a/(2n) 

A - P* a/(4n*) 

Thus, to reduce the aaxiaua velocity while aaiatalnlng the saae 
aooeleratioa a. reduce P by the desired factor and siaultaneously 
reduce A by that factor squared. This analysis Is certainly oaly a 
crude approxlaatloa but should point the way: oae would have to 
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Figure Tether deployed length m • funotlon of tine, 
CoBponentB due to einueold reel*-out oeneuver end active danping 
response shown separately. 
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REEL VELOCITY (m/s) 
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PlguFS 5>10. RmI velocity as a function of time. Ccoponenta 
due to the sinusoid reel-out naneuver and active damping response 
shown separately. 



dtteraint tA« appropriate A for a ohoaan P bp doing a teat oaaa and 
axtrapolatlng linearly aa done above, and poaalbly Iterate on P. 

fe were only able to aake an Initial exploration of tbe releaae 
and daaplng behavior. A nnaber of further poaalbllltlea preaent 
theaaelvea. Aa lapleaented in SCTHOOC. neither the aine-law reel 
aanenver nor the gotlva daaplng ehangea the aaaa need for the 
near-Shnttle tether aegaeat. The error la not large (on the order of 
10%). but ahonld be aooovnted for. 

There la no & priori reaaon that the fixed pre-releaae aanenver 
ahonld be a tlnnaold; thla waa ohoaen to give the de aired groaa 
propertlea with a alaple two-paraaeter faaily. By properly ohooalng 
the aanenver fnnotion (which need not be a alaple analytic 
expreaaion), one no donbt conld rednee the deaanda on the daaplng 
algorltha and on the reel aotor. partlcnlarly In the reel-ont caae. 

The daaplng law alght alao depend on the teaalon tiae derivative 
or Inclnde delay teraa. The good reanlta achieved with the alaple 
law. however, probably do not jnatify thla coapllcation. 

In practice, one would want to Initiate the daaplng law gradually 
to avoid the andden ahoek and reanltlng oacillationa. 

One alght vary the daaplng paraaetera. particularly T^. to follow 

the expected behavior during a aanenver anch aa deployaent. The 

daaper would then prevent excnralona froa the deaired aanenver. One 

could perhapa even create a aanenver by anitably varying T between 

0 

the current and derived atatea and letting the daaper deteralne the 


exact trajectory 
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5.2 Eizlsid Tim til ifi girgplii gihit 

Hara tha triaafar of a payload to a aodarataly high elreolar 
orbit la daaooatratad oslng aa approprlata raal aaaaoTar aad tha 
daaplag algerltha datorlbad la Saetloa 5.1 to atablllta tha 
Sbattla/Tathar/Talaoparator aystaa aftar raleaaa. 

Tha Shattla aaas waa 100 aatrle toaa. tha payload aaaa 9.5 toaa. 
aad tha talaoparator aaaa raaalalag aftar ralaaaa. 0.5 toaa. Tha 

.4 

tathar aaad waa 0.3 ea diaaatar Caviar with a daaalty of 1.5 ga ea 

12 “2 

aad alaatlolty 0.7 z 10 dyaa ca . Tha raqalrad tathar laagth la 
dataralaad by tha ralaaaa ooadltloaa aad la darlvad below. Thla 
laagth la dl ka, with a aaaa of 0.65 toaa. All alaulatloea ware doaa 
with alz aaaa polata: tha Shattla. tha aabaatalllta (talaoparator or 

talaoparator/payload) . aad foar tathar aaotloaa. 

Tha laltlal orbit of tha Shattla waa a 200 to 520 ka altltada 
allipaa (thla orbit eoald ba raaehad with oae bara froa a 200 ka 
elroalar orbit), fa aaaoaad that tha Shattla had baaa lajaeted lato 
thla orbit aad had thaa deployed tha tathar/talaoparator/payload, ao 
that tha orbit eaatar of tha fall ayataa waa la tha praaerlbad orbit. 
Tha prograa which gaaarataa laltlal ooadltloaa for lapat to SCTHOOC 
atarta aa alllptloal orbit at a polat 90** forward froa the perigee, 
la thla oaaa. thla raqalrad aa laltlal altltada of 355 ka for tha 
orbit eaatar. Tha latagratloa atartad at thla polat with a variieally 
deployed tathar aad procaadad aatll apogee, where It waa daalred to 
ralaaaa tha payload. A carta la aaioaat of llbratloa was latrodaoad la 
tha alaalatloa by thla quarter orbit of latagratloa, raaaltlag la a 
■ora raallatle caca. la aa actual, phyaioal aituatloa. tha azlatlag 
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■otloa da* to librati<>a aad ratldaal aotioa froa daplojaeat will htv* 
to b* takoa lato aeooaat ia fiao taaiai tho pr«-r*l«as* aaaoavtr aad 

tlaa of ralaaao. 

Olwaa tba alllptio iaitial orbit for tha orbit oaatar, the laagtb 
of tothar aaat ba dataraiaad ao that whaa ralaatad at apogaa tha 
payload will ba ia aa approxiaataly oiroalar orbit. Aa dafiaad ia tha 
"Stady of Cartaia Laaaehiag Taehaiqaaa Uaiag Loag Orbitiag Tatharu* 
(Piaal Taehaioal Raport for graat NAO-8008. March 1981). tha 
paraaatara ara: 

■ apogaa of pra-ralaaaa orbit (for orbit eaatar) 

Cj " parigaa of pra-ralaaaa orbit 

S ■ radiaa of oiroalar orbit of raloaaad payload 

x^ ■ laagth of tathar froa orbit eaatar to payload. 

P 

The kaowa qaaatitiaa ara: 

- 8898 ka 

- 8S78 ka 

aad froa gaoaatry. S ■ C, ^ x_; 8 aay b* foaad froa aa aqnatioa oa 

* P 

paga 11 of tha raport: 

8* - ♦ CjCj^)/ilCj) giviag 

8 - 8933.48 ka 

* - 8933.48 - 8898 • 33.48 ka. 

P 



Pt|« 94 


ApprotiMtlng th« orbit coator by tho eoator of ■aas» rooalliag tbo 
Mtoot of 100 oad 10 tons, tad aoglootiag tbo aats of tho totbor. tbo 
roqalrod total totbor loagtb io 

I - (110/100) s. - 41.09 ka. 

9 

Tbit oxaaplo atot a loagtb of 41 k« oxaotlj. Tbit it tbo ttrotehod 
loagtb of tbo totbor ia oqailibriaa at tbo iaitial poiat oa tho orbit 
and it a roqairod for iapat to tbo iaitial ooaditioa roatiao. It 
woald baro booe aoro prooito to adjatt tbo atrotebod length of 41 ka 
at apogoo back to tho loagtb ia tho gravity gradioat at tbo iaitial 
poiat. bat tbo difforoaeo ia taall ooaptrod to othor aogloctod offoott 
taob at libratioa. 

Tbo dotirod olliptieal pro-roloato orbit aad tho totbor loagtb 
vbioh will allow roloato iato a oirealar orbit bavo aow booa 
oaloalatod. Tho fiaal ooapoaoat aoodod it tbo rool aaaoawor aad 
tiaiag of tbo roloato. Booaato of diffioaltioa ia roaowttiag tbo 
prograa to goaorato olliptieal iaitial ooaditioat. tbo rool aaaoavor 
wat iavoatigatod toparatoly. atartiag with aa iaitial eiroalar orbit, 
aad lator eoablaod with tho iaitial ollipto. Tbit it aa aoooptablo 
proeodaro for a firtt approxiaatioa. bat earo aatt bo takoa to oporato 
with aa orbit of tbo taao boigbt at tbo iatoadod roloato: Tho totbor 

toatloa vtriot tt tbo gravity gradioat* vbleh variot at Tbo 

toaaioa la a totbor at 910 ka altitado that diffora froa that at 200 
ka by (4891/4971)*'^ • 0.S47. l.o. 13% lott. Otbor offoott. 
otpoeiLlly libratioa. will offoet tbo toaaioa. bat tboir eoatribatioa 
it aiaor (tbo total dlfforoaoo aotod whoa tbo ollipti. oaao vat raa ia 
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f«il was 0.8S6). Sinoa tha raal nanaavar la datignad to rednca tha 
taaalon by toaa 95^. a aanaaver snltabla for 200 ka altltuda will ba 
■oat oaaaitabla for aaa at 520 ka - Indaad, it will oaoaa tha tathar 
to |0 alaek. With this eaatioa in aiad. tha followiag paraaatara wa 
foaad to ba appropriata: 

Raal Naaaavar: A aia(2n t / P d) 

A “ -182.24 aatara 
P " ISO saeoada 

6 • - n/2 for gradaal iaitiatioa of aaaaavar 

stoptiaa ■ P/2 for gradaal taraiaatioa 

Aotiwa Oaapiag Algoritha: A(t) > RlKtl-T 1 -8 A(t) 

0 

■ 0.16 z 10^ dyaet 

^ ■ 0.14 X 10 ^ ea dyaa ^ aso ^ 


8 - 0.047 aac"' 

Tha followiag foraalaa. baaad oa aoaliag argoaaata with coeffieiaata 
dsrivad froa a aaapla oaaa, provida good atartiag poiata for tha 


acth 


8 ■ 0.041 aae ^ [d/0.3 ea1[x/80 ka1 [a^/0.5 toal 

K ■ 0.16 X 10 ^ ea dyaaa ^ aae ' [d/0.3 oa1 ^ [i/80 ka]'^^ 


[a^/O.S toa] 


- 1/2 


T • 0.:86 X 10 dyaaa [x/80 ka] 
0 






924 toa J 


0.848 toa 



althoagh tha optiaoa taaaioa to aaa ia baat derived by aeaaariag 
tha aotaal taaaioa in tha daapad. ralaaasd atata ia a aort of 
bootatrap proeadara. Wars, d ia tha tathar diaaatar, s ia laagth. a^ 
tha aabaatallita aata. aad a^ tha tather ataa. Tha rslaaaa tiaa waa 
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ohoa.-n ai tha moaaat when tha tantios at tha sobsatallita (payload * 

talaoparator) was tha aaaa aa tha oltimata aqnillbriam teaalon at the 

talaoparator after ralaasa. Aa with tha ta&aioa T at the Shuttle, a 

o 

flret approilaatioa to thla desired taaslon is 

’^subset " dynes (a/80 kal (b^/ 0.5 tonl 

whioh ainst than be iterated to nse the actnal relaxed tension after a 
sneoessfnl, but sub-optiaal, release is achieved. The release tiaie 
derived by this siethod is 55 seconds. One final point to note: this 

man,»aver uses a strict reel-ont control law, without an initial 
reei-in. This leaves the released payload with a slight outward 
radial velocity relative to the orbit center; thus, to get a circular 
orbit (zero radial velocity), the release must occur slightly after 
apogee, when the system has a net inward velocity. 

These components can now be "patched” together. First, an 
elliptical orbit was generated beginning with a vertically deployed 
tether 90*’ forward from perigee. This orbit was followed until 
somewhat after apogee (1/4 orbit) with no release. A sample run of 
the reel maneuver without release was also prepared. Programs were 
written to read the SCYHOOK simulator output and print tables of 
various quantities at each output interval, specifically, the radii 
and velocities in the radial and in-plane directions of the Shuttle 
and subsatellite, the circular velocity for these radii, and the 
velocity and position of the subsatellite as perceived in the 
coordinate system centered on the Shuttle and with axes in the radial 
and in-plane directions (all out-of-plane components were zero). From 
these tables, the time at whioh the system's inward velocity would 
balance the release velocity relative to the orbit center was 


- "3U 


r 


•stlMt«d. This it tbt first s 9 proxiauitioa to tht rtlettt tiat. tad 
by sabtrsotiag SS ttooadt from tbit tiat t tttrt tiat for tbt rttl 
aaataTtr vat obtaiatd. 

Soat dttailt, taeb tt tbt tfftett of libratioa, atrt atgltottd ia 
dtriviag tbttt start aad rtltttt tiats. Ibt atzt stop was to rtfiat 
tbt tttiaatts. To redact ooapatation tiat a SCTBOOC iapat filt vat 
ortattd. oa tbt basis of tbt SCfHOOK oatpat froa tbt iaitial 
tlliptloal orbit, vbieb allowed the start of tbe iategratioa shortly 
before tbe aaaeaver start tiae derived above. Soae leeway was left so 
that this start tiae eoald be varied. Several orbits were raa with 
differeat aaneaver start tiaes and no payload release. For each start 
tiae. a release tiae was defined at that tiae vbea tbe total radial 
velocity of tbt sabsattllitt (payload) was stro. Tbtn, by tabnlatiag 
tbt reltast tiat rtlativt to tbt start tiat (this thonld be 55 seconds 
for optiaal post-release tetber behavior) aad tbe difference between 
the actaal la-plaae velocity aad tbe circalar velocity at the 
sabsatellite radina (which shoald be aero for a circalar release) at a 
fanctioa of tiae. it was discovered that tbe first qaaatity varied 
ancb aore rapidly than tbe others, so tbe relative release tiae was 
obotta to deteraiae the start tiae. aad beace tbe tiaiag of the entire 
aaneaver. 

That, by beginning with a deployed tetber in elliptical orbit, 
ooadaetiag a tensioa-rednciag aaneaver aad releasing tbe payload at a 
tiae carefally cbosaa to control tbe tether's post-release 
oscillatioas and to leave the payload with aero radial velocity, tUe 
payload orbit was placed ia an alaost perfectly circalar orbit. The 
final eeeentricity was 1.7 a 10**^ which was deeaed scceptable. 
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No attoapt was aada to aaalyaa tba daptndaaoa of tlia flaal orbit 
oa tba otbor paraaatara available, aor to aodify thoae paraaatera to 
aobiOTO a aor# aaarly eiroalar orbit. The priaary paraaatera 
available are the laitial orbit paraaetera (aay the aiaiaaa aad 
aaiiaaa heighta) aad the tether leagth. The teleaae aaaeaver haa 
relatively little iaipaot. the velooitiea lavolved beiag aaall (6 
aeter/aee) ooapared to the radial velocity of the elliptical orbit 
evea a aiaate froa apogee (12 a/a). How aight theae paraaetera be 
adjaated to prodace a eiroalar orbit if the theoretically derived oaea 
had aot beea adeqnate? Froa a coaputatioeal ataadpolat. the aioiplect 
way voald be to adjaat the apogee height, redaciag it to laoreaae the 
payload releaae velocity relative to eiroalar velocity. The release 
aaaeaver woald aot aeed to be reooa^ated. oaly ceded. la a phyaioal 
aitaatioa. the orbit ia aore likely to be givea aad thea oaly the 
tether leagth eaa be varied. Thia voald reqaire recoapatiag both the 
iaitial orbit aad the releaae aaaeaver. ladeed. aa ideal ooatrol 
algoritha for releaae woald alao coataia a relatively brief (perhape 5 
to 10 aiaatea) deployaeat/retrieval phaae for aiaor adjactaeata to the 
tether leagth before atartiag the teacioa redaciag pre-releaae 
aaaeaver. To achieve eiroalar releaae at a preaoribed height, two of 
the three paraaetera (apogee height, perigee height, tether leagth) 
aaat be deteraiaed. leaviag oaly oae free paraaeter to aeet other 
aiaaioa reqoireaeata. 

Oae each aiaaioa reqaireaeat voald be for the fiaal orbit of the 
ayatea to be high eaoagh after releaae to avoid re>eatry. SKYHOOK waa 
allowed to iategrate the ayatea after payload releaae well pact 
perigee. The Shattle'a loweat altitade waa deteraiaed to be IdS ka. 
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dangtrously low la a praotloal sitnatloa. Thara was no ratriaval of 
tba tathar, bat thia wonld calsa the orbit only alightly. To avoid 
thla aitnatloa tha initial orbit wonld aaad a hlghar parlgaa, 
raqalrlag two ahattla burns for injaotlon. 


To Sanaa rite: 

Shuttle: 

100 metric tons 

Sab satellite: 

O.S-ton payload 

0.5'ton teleoparator (remains after release) 

Tathar: 

Length: 61 km 
Mass: 0.6S ton 
0.3-om diameter Kevlar 
Initial orbit: 

Perigee altitude: 199 km 

Apogee altitude: S19 km 

Boeantrielty: 0.0237 

Started 90*^ after perigee, with tether fully deployed 

Reel maneuver begins at 1386 seconds after start: 
parameters 

Aaiplitude: 182 meters 

Period: ISO seconds 

Half-period sinusoid maneuver, shifted so that onset 
and stopping are gradual (a cosinusoid). 

Pure reel-out maneuver. 

Release payload at 1441 seconds. 

Active damping of tether begun at release. 

Final orbit of payload: 

Eccentricity: 0.00017 

Altitude: 57S km 

Pinal orbit of Shuttle: 

Perigee altitude: 165 km 

Apogee altitude: 514 km 


The tension and a side-view of the system are plotted for t * 0 


to 2000 seconds in Figures 5-11 and 5-12. The tension behavior near 
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Figure 5-11. Tether tension froa the initially deployed 
configuration (t - 0) through the release and daapii« (t -1441) until 
t “ 2000 seconds. The detailed aission profile is described in the 
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tb« tia« of rtlotto is plottsd in Figars S-*13. with sn sxpsaded tlao 

soils. Figure S-14 gives s soheastlo represeatstlon of the psylosd 

orbital trsasfer. 

S.3 Psylosd Aoqaisitioa by s Tether Deployed froa s Spsoe Ststloa 

The tether behsvior whea s psylosd is saddealy soqaired by s 
tether deployed dowawstd froa s Spsoe Ststioa ia aoderstely high orbit 
was iavestigsted. 

For this speoifio ease the Spsoe Ststioa was plsoed ia s 500 ha 
sltltade oiroalsr orbit. The psylosd was iajected froa so iaitisl 200 
ka psrkiog orbit, giviog a 200 ka perigee for the psylosd delivery 
orbit. The tether leogth sod the apogee height of the psylosd orbit 
were theo sdjasted so that 1) at apogee, the psylosd is st the ssae 
height as the teleoperator eod of the tether sod 2) the velocity of 
the psylosd ia orbit st apogee is the ssae as the tethered 
teleoperstor velooity. 

If R is the ststioa* s orbital radios sod L is the tether leagth 
thea the teleoperstor velooity is 
V* - [(R-U/Rl* OM/R 
or 

- (1-k)* (OM/R) 

wbtro X ■ L/R ii wbat wo want to finds 

Tho payload orbit's porigoo is spooifiode The saall paraaoter 8 is 

introdneod so that the perigee radius is (l-*8)Re The apogee is given 

by condition (1) as (l-X)R and hence the seaiaajor axis is 

• • tl*(X^8) /21Re Theng epplying • OH [2/r - X/al at apogeci r ■ 

(l-X)R. the psylosd velooity is V * ■ (QM/R) [2/(l-X)-l/(l-(X+8)/2)] .. 

F 
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PAYLOAD TRANSFER TO CIRCULAR ORBIT 



Plgur* S-14. A 8ob«Mtlo view of the verioue orblte dlscueeed in 
the text. The beevily drevn portione repreeent eotuel trejeotorlee. 


original page is los 

OF POOR QUALITY 

Coaditlen (2) raqniraa •qutUg V_ tad V_, Uadlai aft«r tea* 

P * 

■aaipalatloa to 

tl-(X+d)/2l - 1-5 

Tkit it aa tqattloa to 5t tolvod fot X. tht totltd totbox Itagth, 
glTta tht partatttr 6. tht totltd difftrtaot httvota tht Sptot 
Stttloa tad ptplotd prt-iajtctioa orbital radii. Siaet X tad 5 art 
tatll. t firtt-ordor tolatioa girtt X«5/7. Tht aqaatioa «tt tolvtd 
aaatriotlly for 5 i t0.11 tad tht ratio (7X/5) to tht firtt-ordtr 
tolatioa it plottad ia Figaro 5-15. For thit otto 5 ■ 

(500-200) /«378’f500) - 0.0435. which it iatorpoltto to (7X/5) - 1.0275 
giwiag L ■ 44.035 ha. Thit it tht loagth froa tht orbit ctattr to tht 
toltoptrator. Uaiag tht atttta girta bolow tad t etattr-of-aatt 
approxiaatioa. t total tathtr Itagth of 44.07 ha ia obttiatd. 

Ia taaatry. tht tyttta paraatttrt art: Sptot Statioa orbit. 500 

ha altitado tiroaltr; payload orbita. 200 ha ptrigtt tad 455 ha 
tpogtt: Sptot Statioa attt. 1000 attrio toat:Ttthar paraattart. 44.07 

ha loag. 0.3-oa diaatttr Ctvlar:Ttltoptrttor attt. 0.5 toa:Paylotd 
attt. 9.5 toa. 

Figart 5-15 thowa tht ttaaioa rttpoatt whta tht payload it 

taddtaly toqairtd at t * 10 ttooadt. Nota tht ttroag ttaaioa 

oteillatioat with ptriod tboat 185 ttooadt. with atziaaa ttaaioa aboat 

3.2 z 10^ dyata. By tpplyiag tht activt dtapiag algoritha diteatttd 

ia Statioa 5.1 (abort) tht ttaaioa rttpoatt of Figart 5-17 it 

tehitvtd. Tht dtaptr wat taratd oa at tht aoBtat of aoqaititioa. with 

paraatttrt C ■ 2.55 z lO”^ ea dyatt”^ ttc”^. 8 ■ 0.012 tto”^ tad T ■ 

0 

1.5 z 10* dyatt (thotta tlightly Itat thaa tht tzpttttd 
paat-atqaititioa tqailibriaa ttaaioa of 1.58 z 10* dyatt to rtdaot tht 



ORIGfNAL PAGE It 
OF POOR QUAimr 


Figure S"^1S. X ve. 6. X is conputed froB the equation 
[l-X]^U-(X*^)/2]-l^ for the paraoeter 6 with range 0 to 1. X is 
the sealed tether length (L/K) and 6 is the ainilarly scaled 
difference in height between the station orbit and the payload parking 
orbit. The ratio of tk^ aotual solution X to the first- order 
approxlaation 6/7 is plotted. 
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first tension msiimun). Tlie system reaches a state of virtual 
equilibrium in about six minutes and the maximum tension reached, 

g 

about 2.15 X 10 dynes, is only about 28% above final equilibrium and 
67% of the maximum in the undamped case. 

For future comparisons, the damped case of Figure 5-17 is plotted 
to a different scale in Figure 5-18. Figure 5-18(a) shows the tension 
vs. time w!iile Figure 5-18(b) is a side-view plot (radial vs. in-plane 
components) of the tether state at two-second intervals, with the 
spacing between mass points reduced to 1/2 km to emphasize the 
motions. Without this scale expansion. Figure 5-18(b) would appear as 
a motionless '^picket fence**. 

It would be desirable to reduce the tension overshoot from that 
found for the simple damping-at-acquisition case. As a first 
approach, the damping algorithm was turned on before acquisition, 
allowing it to increase the tension somewhat. The tension vs. time 
plot for such a run, with the damper initialized t - 0 and the 
payload acquired at t == 9 seconds, is shown in Figure 5-19. The 

damping behavior is substantially improved and the maximum tension is 

8 8 
only 1.95 z 10 dynes. 16% above tbe eqoilibrinm of 1.68 z 10 dynes. 

g 

compared to 2.15 z 10 dynes for the simple damping case. 
Unfortunstely, this method is not very stable: Figure 5-20 shows a 

case in which acquisition did not occur until 20 seconds after the 
damper was initialized at t » 5 sec. The tension plummets after about 
10 seconda of damping and the tether goes slack soon after 
acquisition. These two ezamples show that the peak tension, hence the 
tether strength needed, can be reduced by suitable pre-acquisition 
maneuver, but that slight changes in the manenver can have drastic 
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Figure S-19. Post-aoqulsltlon tension variations are reduced by 
turning on the active damper at t ■> 0 and delaying acquisition until t 
■ 9 seconds. 
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MASS SYM30L 



Flf. 'o 5*20. Tanalon va. time when tho payload Is not acquired 
until 20 (teconda after tho daunper la started (at t •> 5 seconds). The 
tension had already begun to decrease when acquisition occurred and 
went alack at t •• 31 aeoonda. (Note the expanded tine scale In this 
pxot.) 
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offocta. Rxperlmeatatton with a fixed sinusoid maneuver, as used In 
the previous release etudles. would be worthwhile. This situation 
would probably be less sensitive than the use of the active damping 
algorithm to generate the maneuver. 


The payload transfer orbit la likely to be slightly In error, so 
that the tether will suffer additional disturbance due either to 
maneuvering of the teleoperator or to acquisition of the payload with 
a non-zero velocity. Several cases were run in which the payload had 
a small residual velocity (3 m sec at acquisition and with the 
velv^clty vector in various directions. The tension and side-view 
plots are given in Figure S-21 (velocity In-plane) and Figure 5-22 
(velocity out-of-plane). We also show the front-view (l.e. radial vs. 
out-of-plsne components) In Figure 5-23 (velocity radial, directed 
away from the station) and Figure S-24 (velocity radial, directed 
toward the station). In all of these cases the payload was acquired 
and the damper Initiated at t ■ 0. Also, the radial component scales 
in all oases are expanded by a factor of about 18 to show the motion 
clearly. 

The fundamental lesson of these cases Is that the situation is 
stable and the damping algorithm is able to cope adequately with these 
disturbances as well as the simple acquisition with zero relative 
velocity. 


The tether response la all cases is much as expected. The 
tension response la the In-plane and out-of-plane oases. Figures 
S-21(a) and 5-22(a), is very close to that in the zero-velocity cose. 
Figure S-18(a). which is understandable since tension is force along 
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Figure 5-21. Plots, for the case in Figure 5-18, but with 
acquisition occurring with a 3 ateter/sec relative velocity in the 
in-plane direction along the orbit, (a) Tension vs. tine, (b) Radial 
vs. in-plane behavior. 
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Figure 5-22. As for Figure 5-21. excej-t that the velocity is in 
the out-of-orblt plane direction, (a) Tensi .’n vs. tine, (b) Radial 
vs. in-plane behavior. <c) Radial vs. out-of -plane behavior. 
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Figure 5-23. As for Figure 5-21. except that the velocity Is 
radial, directed away froo the station (downward), (a) Tension vs 
tiwe. note the tension variation is greeter than the other cases. 
Radial vs. in-plane behavior. 
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Figure 5-24. Is for Figure 5-21, except that the velocity is 
radial, directed toward the station (upwf’d). (a) Tension vs. tiae 
Mote that the tether segaent nearest the payload goes briefly clack 
This aay be an artifact of the luaped-oass aodel. (b) Radial vs. 
in-plane behavior, the dashed lines ir llcate slack tether segments. 



Mfi5S STNSO^ 




no'n i; 


0<M-XO^ 
— rvi m a* i/w£: 


OmGfNAL PAGE IS 
OF POOR QUALITY 




fNl 

I 

LH 

a> 

u 

3 

OG 

U, 





Page 124 


thf tether. Tn these two oetes» sabstential swinging of the tether is 
seen. Figores S-lKb) end Figure 5-22(b), in the direction expected. 
When the welooity ii rsdinl. little swinging of the tether is seen. 
Figures S-23(b) end Figure 5-24(b)« but the tension response is 
sabstsntiil ly different fron the zero velocity esse. With the 
velocity iwey-from-stst ion. the Initisl bounce is substsntislly higher 

g 

(2.S5 X 10 dynes, compered to 2.15); while with velocity toward 

g 

station, the tension peak is somewhat lower (2.05 x 10 dynes). This 
substantial increase in tension with • 3 m sec ^ velocity 
awsy-fron*station indicates that this direction is clearly less 
desirable. 

*rho voloolty-toward-station case (Figure 5-24) shows a brief 
initial period of about 6 seconds when the tether segment nearest the 
payload goes slack. This is rather surprising since tension 
disturbances propagate through the tether very rapidly compared to the 
3 m sec ^ acquisition velocity (see Figure 5-20) in which the final 
segment, at 40 km distance, shows substantial response a mere 5 
seconds after the damper is initiated. The slack response is probably 
an artlfsct of the lumped mass model used. This should be explored by 
ohsnglng the number of mass points to see if it is only the outermost 
segment which goes slack in all cases. The impact of a 
constant-velocity ram on a one-dimensional tether could also be 
inveetigated theoretically. 

One countcr-intnit ive response of the tether is illnstrated in 
Figores S-18(b). Figure S-23(b). and Figure 5-24(b) for the 
zero-veloeityi •way--fro»-statioa end toward-stet ion cases, 
respectively. When the extra maas is attached to the tether. 
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Iccreteing the tenelon, the expected response, et least initially, is 
for the tether length to increase. Instead, it decreases (but note 
that the aasa oloseet to the station initially moves farther from it). 
This is even more pronounced in Pignre 5'*23(b) where the velocity was 
directed so as to stretch the tether. While when the velocity et 
aoqnisitlon should have compressed the tether (Figure S-23(b), instead 
It responded by lengthening. This probably represents an Interaction 
with the active damping mechanism and requires further study to be 
understood more clearly. 

To summarize, these studies have shown that payload acquisition 
by a tether deployed downward from a Space Station in high orbit is 
feasible, that the damping algorithm described in previous reports 
suffloes to reduce the initial tension increase and eliminate tension 
oscillations within a few minutes, and that the procedure is stable 
against acquisition with non-zero relative velocity. But further 
investigation is Indicated to develop pre-acquisition maneuvers to 
reduce the tension bounce, to resolve the question of whether the 
tether actually does go slack if the payload is acquired with a 
vfloclty toward the station, and to understand some of the puzzling 
system responses on a clear and intuitive level. 

Qr.feiUl >I.'iaP.L!ll 

5.4.1 Introduction 

In this section we consider the general problem of orbit change 
and control by varying the shape of a spacecraft by internal forces. 
This docs not mean that external forces are not operating. The 
external field remains Newtonian but the gravity gradient forces 


Ir 
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ooablaed with the shape variation of the spaoeeraft oanse a variation 
of the total force acting on the center>of-maas and, therefore, ^anse 
a variation in the orbit. Since the external forces are central, the 
total angnla momentna does not change. What may change is the 
orbital energy, which may either increase or decrease. 

S.4.2 Tether Model for Orbital Pumping Stndy 

Here la a very simple mechanical system for modelling the 
mechanism we intend to investigate. 

Consider two eqnal maeses m at the end of a tether in orbit: a 

dumbbell (see Figure S*-25). The analytical model of the tether is 
assumed to have the characteristics of a bar capable of withstanding 
compression as well as tension. In reality, we are interested only in 
the tethered "dumbbell*' configuration and. therefore, only in one 
class of its possible movements- those for which the tether is always 
in tension. If the orbit is equatorial (about an asymmetric primary 
body) or if the primary is spherioal, the external forces on the 
dumbbell are central and consequently the total angular momentum of 
the spacecraft is preserved. If a and e are the somimajor axis and 
ocoentrioity of the orbit of the center of mass, wc have that in any 
ovolntionary motion which does not effect the orbital angular momentum 

a(l-e^) - a* (1) 

If by internal control we change the orbital energy, which is to say, 
the semimajor axis a, we also change e. If a increases, e also 
increases, and vice versa. From 1) it is clear that a* is the minimum 
value of the semimajor axis of the orbit and corresponds to a circular 






Figure 5*25. Co*ordlnate relationship in dumbbell configuration 
used in orbital pumping study. 
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orbit. Nitvrtlly. a* la tba taallatna raotiia of any orbit. 


Hanriatloally. whaa tba dnabball ia la orbit, fraa aad foroad 
libratloaa will oooar. fa aaauaia that tba taaaioa ia tba tathar nawar 
goaa to saro. Aa an azaapla, lat aa aaanaa that the orbit baa amall 
aooantriolty and that tba libration aaplitadaa are aaiall. The tathar 
ia alatoat alwaya allgaad with the radial dlractioa. When the ayataa 
ia at parigaa, tba tanalon ia the tathar ia roughly given by: 


Ra^za 

% ■ 

If 2z is tba length of tba tathar. whan the ayataa la at apogee the 
tanalon ia 


( 2 ) 




zm 


M V «ani 


(» 


If wa let a aegaent Az of tether be pulled out at perigee by the 
gravity gradient tanalon againat the torque of an eleetrlo generator 
and wa pull in with the aaaa aotor at apogee exploiting the generator 
power for raeycling tba aaaa tathar aegaent. wa will have tranaforaad 
orbital energy into alaotrio energy (for axaaiple) or into thamal 
energy approziaataly equal to 

lz(t^ - x) • 18 ag^ -T- za Az 

P t ^ 

which mty bs writtoo 


(4) 


Az (t - T ) - 18 a li a (5) 

pa • a* 

Since the orbital energy la |i/2a. whan the eccentricity is pasalng 
froB 0 to e. while the angular aoBantua la praaarvad. the variation of 

energy is 
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» _ lJ2 + -Jf g . 


AE - 

• 2.(l-e*) 


• 1-.* 


““ t : 
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( 6 ) 


Thetafora, for amtll aoeantrtoitlaa. tha fractional variation par 
orbit ii obtainad by dividing 5) by d) . fa gat 

AE/orbtt ^ jg tAa(l-a^) 


AB 


aa 


(7) 


In partionlar for a ■ 100 kai, Az ■ 20 km, a ■ 7000, a = 0.03. 


. J.4 . I0-* 

< 

Notloa that thia laat valna eorraspondt to tha valua 0.03 of a. 
Tbarafora, tinea tba fractional variation of energy par orbit it 
invaraaly proportional to a tha procaat la not uniform. 


( 8 ) 


In fact, «a may vrita, from 1) 

E - ns . a 

a a* 

and bacanaa of 5 ) 

' • -2aAa ^/orbit ■ 18 a zAz/orbit 

orbit a* ^,2 

and, tbarafora. 


(9) 


( 10 ) 


Aa/orbit - 18 

a** 


( 11 ) 


While other perturbat iont may give abort periodic variatiuna, tha 
atmoapharic drag ia dacraaaing tha apaoacraft orbital energy and 
angnlar momantum. Aa a flrat approximation, for amall aocantricltlaa 
tha variation par orbit of tha energy la given by 


3l. 

orbit 


2« a Ap(a)v^(a) 


( 12 ) 


while tha variation of tha angnlar momantnm la given by 
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0 1 

AL/orbit ■ * Ap(a)v (a)T 


Sine* for a olroaUr orbit 

^ a^a' ■ »». T 


2 3 


2 ^ .2 3 _ *2 . 4n‘a 


ire have 


(13) 


(14) 


AE^/orblt • C^nttAp(a) ■ 2 CQf(Ap(a) a 

AL/ orbit ■ nCj^A p(a) a pa (14*) 

from A), and 12*) wa have 

^^n/ofblt 2CDnAp(a)a 

< ‘ ,,2 (IS) 

if «* aaauai* A • 10*^-10^ oa^, C^j ■ 2.2. ■ ^ 100 z 10^ graai*. a ■ 7000 
I 10^ oa^, p ■ 10 grami/offi^. * ■ 0.01, we have 

( 15 ') 

4B° 

e 

to be eoatparad with i ) . 


Similarly 

• 5 . 10-‘ (1« 

L 

f'or a dumbbell ayatem flying inalde the range of 400 to 1000 km, the 
almoapherlo drag ia probably negligible for a realistic configuration 
of the syttom. The preaent technique in this situation may therefore 
be useful in controlling the orbital elements a and e in the preaeuce 
of small atmospheric drag and any other small perturbations such as 
radiation presanr*. 
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S.4.3 Analjrtloal Model tot Nuatrloal Slanlatlon 

A nuaerioal tlaolatloa will give all th* *l*a*iit* for a praoia* 
•valvatiofi of th* taohalqa*. la th* followlag paragraphs w* ahall 
ooasldar a r*latlv*ly alapllflad analytloal aod*l la order to be able 
to perform a prellalaary ayatea aaalyala which, however, will 
eveataally reqair* oaaerleal lategratloa la th* aoat latereatlag oaaea 
(large eoeeatrlolty) . 

For eatabllehlag the dlffereatlal eqaatloaa of aotlon of the 
ayatea of Figaro S^lS. w* latrodno* the followlag Lagraaglan 
ooordlaatea 0. r. A (ae* Figaro S-2S). If w* also ooaslder aa aa 
ladepeadeat variable th* leagth of the tether 2z we may then ooaslder 
th* foar degree of freedom system with ooordlaatea 9. r. A. s. 

W* first flad th* klaetie energy 

1 1 /4*V\ 

T - J 2 " ^*2 ^ 2 ^ 

la Cartesian coordinates and la Lagranglaa coordinates 

t - m (f* ♦ r*4* ♦ (18) 

Th* force fanotloa of the central field of gravitational force 
centered la 0 la 

u . liS! + m (19, 

* (t^+i*-2rsco*A) (r*+a*e2racoaA) 

If th* tether Is vlsooelastlo w* have two more forclag faactloae. the 

elastlo force fnaotloa: 

®c • - f 


( 20 ) 



and th« dlaslpattoa fnaetioa 
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R - -2 


Di 


2 


( 21 ) 


Oth*r for*** ooald ba appliad to tha ayataa aach a* aarodyaanlo 
foreas. radlatloa pcataara. and potaibly alaotrodyaamlo foroat: 
ho«*v*r< «* will take la ooaaldaratloa la tbls prallmlaary nodal only 
tba foroaa «* faav* Hated above. Natarally tbe nodel laolndea tbe 
oaa* of ooatrollabla laagth da of tbe cable. Thla coatrol will Inply 
lateraal forcaa, and related work, and exchange of energy between 
orbital and rotational energy and Internal energy of the ayataai. 


Control forcea nay be alao laclnded explicitly. la partlonlar 
two aqoal and oppoalta forcea nay be applied to and P 2 . Tbe 
oorraapoadlag work nay b* written 

6 • 2fdz 2Tz6(e-«-d) (22) 

whar* f la the longitudinal and t the tranaveraal conponant. If t « 0 
there la no variation of angular aioaiantun of the ayatea. Limiting 
onraalf to tha force* explicitly ooaalderad we may writ* tha 
Lagranglan equation of not loo 

2 2 

2n(r-r«*) - • -ipi ♦ 3|i5*_ - 9„5L- coa^d 

• r r r r 

2n ~(r^O-*-z^(0+d) 1 » 2tz 
dt 

2 

2n l*^(R+i)l • -6 I ”- aladeoad + 2 zt 
dt 3 

r 

2nt*-z(«+i)*l - -2k(*-z^)-20i -2u5i(l-3coa*d)+2f 

0 3 




( 23 ) 
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W« ooA«id«r la tha praiaac parairaph aoma apaolflc oaaat. Plrat 

wa oontldar tha alaplaat oaaa x *0, k * • and tbarafora 

a * t ooaataat. Tha apataa 23) baooaaa 
0 


J.,9* . JL ♦ ili - ifi 


* r* 


+ i*(e+i) - L 


2 •* af 9 )1Z 

t*(e+2) - - aiald (23' ) 

0 3 

r 

Thla diffarantlal ayataai aodala tha drag-fraa rigid dumbbell. In 

partloalar. tha ooapllag batwaan llbratlon of tha dombball and orbital 

2 2 

action of tha oaatar of gravity appear axpllcltly. If z << r^ wa 

0 o 

any decouple tha ayataa and nalng perturbation tachniqnaa atudy la 
flrat approzlaatloa how tha action of tha duabball abont Ita 
cantar^of-gravlty affaota tha action of o.g. and vice veraa. In the 
aaoottd oaaa. wa control 6 to tha value 6 " 0, applying two oppoalta 
foreaa. given by 


T • (z5 + 2zO)a 


(24) 


It la evident la thla eaaa that tha total angular aoaantua la not 
praaarvad (aae aquation 23.2). fa treat thla oaaa bora Incidentally; 
but it la natural that even though tha total aztarnal force la taro, 
the aztarnal torqna la not. 

Tha force 24) la aaall only if i and 0. 0 arc aaall. Now for a 

a *3 ** *2 

low orbit tha order of 0 10 rad/aao. the order of 0 la 2a0 . 

“2 

Aaaualng a “ 10~20 a/aac and a ■ 10 tha acealaratlon la of the order 

"3 -3 

of 10 g. which aaana that t aay reach value of 10 ag. a quite large 

force la particular If tha ayataa haa to be operated for long perlode 
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of tl«« while tho length of the tether t ie controlled. 


In thle ctee the equation* of notion becoae 


2‘ 

re - L 


f - t •* - - i- 

T t 

; - .4’ . - i ) - D i . 4 Hi .f 

■ Ob 3 

r 


d*l/r 


( 25 ) 


( 26 ) 


end following the Blnet fomnle. 

SQ^ L ^ L r 

ov o o 

Sinot z it i pontrol funotion (of 0 ) tht problen may bo aolvad only 
nuBtrioally la tbt oaat of larga eeoentriolty of the rtftrenoa orbit* 
Por tmall aooontrlolty on# nay nta a perturbation teohniqne* We 

attame 


•(!-* ) 


( 27 ) 


and 26 ) b*ooai*« 


de* 


•*■ » - 


« 3 2 

[i+, 

ooa(e-u 1 

L 6 
o 


If we aaanai* 


( 28 ) 


a ■ * (1+0 ooa(e-e*)l 

0 


Equation 29 ) beooaiea 


- 2,12 

X • ^ ~ [l-«-aeot(e-e*) ] [l-t-ecoaie-w) 1 

de* K'" 


( 29 ) 


( 30 ) 
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•nd 

d*. 

— - + i»k[l*-<»ooi(0-®*)+*oo»(0-ii) 

« ^oa(t« -«• -a)) Ol) 

2 2 

A ••ntrtl solution of this •qnatlon contains (a) a constant tera 

*- ■ kt“ oos(®*“»)+11 (31') 

1 2 

(b) seonlar taras 

I, - ^ e sln(O-0*), X - k e sln(9-») (31" ) 

2 2 3 2 

(e) hl|htr fraqntnoy term 

- -k ~ cos(29-9«-w) (31'" ) 

Ths oorrtotlon on 0 should ba dona a postarlorl oonsldarlii 25). 

• L tba constant a. 9*. s ara part of tha control 
0 0 

function. 

Tha third aorc natural and raallstlc case is to have t» 0. f*0. 

In this ease, tha aquations baooaa aora coaplaz and tha aodal 
rasaablaa a satallita In orbit affactad by tldas raised by ths prlaary 
on tha duabbsll. Tha dynaaloal aquations ara practically 23) 
where T*0 and f*0. For anhanclni tha tidal effects, the elastic 
constant k Is aada vary saall and a critical das^ini constant Is 
salactad. 

Nora eoaplax Is tha fourth case where Isngth changas are controlled by 
raallug operations la addition to being due to viscoelastic behavior. 
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For all foar aodala aoaaldartd abovt. oaljr auaarloal iatagratloa 
will glva lataraatiag raaalts la tha aaaa of high aooaatrielty whloh 
appaara to poaaibly ororlda aoaa aaw taehalqoaa for aodifyiag tha 
orbital alaaaata of tha e.g. by lataraal foroaa. For aaall 
aooaatrielty orblta. whaa tha ataoapharle drag la aagllglbla. tha 
third oata t"0. f*0 haa baaa daalt with la a noabar of papara. tha 
■oat laportaat ona la tha papar by P. Ooldraleh. "On tha Booaatrlolty 
of Satalllta Orblta la tha Solar Syataa" (O.R. of Royal Aatronoay 
Soolaty. vol. 126. 1963. p. 357). la that papar It la ahown that for 
a Boa-aplaalng aa tall It a. la oar oaaa for a llbratlag iuabball. tha 
tlda ralaad oa tha aatalllta aay ovarooaa tha affaet of tha tlda 
ralaad by tha aatalllta oa tha prlaary. Ihlla tha foraar taada to 
daoraaaa tha aeeaatrielty. tha lattar taad to iaoraaaa It. For a 
larga daabball. avaa oaa of tha alaa of a larga apaea atatloa. whaa 
tha dlaalpation factor la aada larga aaoagh. tha acoantrlolty will 
eartalnly daoraaaa If ooatrol ayatoa la oparatod on tha daabball. 
Tha aquation of motion will ba aquation 23) whore wo put x t * 0 
and lat tha orbit ovolva undar aatural foroaa. If tha dlaalpatlva 
tide ralaad by tha Barth oa tha daabball la larga (amall rigidity and 
low quality factor of tha tathar. or tha laaartlon of a aprlng-daahpot 
ayataa with proparly daalgaad oharaeterlatloa) wa will gat 
olroalarlaatloB of tha orbit. 

For aora elarlfloatlon. wa will eoaaldor raaonant apla*orbltal 
ooupllag dlffaraat froa ona-*to**oaa ratio (whloh aaaaa a aotloa 
dlffaraat froa tha daabball aahlag oaa rotation aud llbratloa during 
ona ravolution). Tha eonstaat aagular aoaaatua la ralatad to tha aua 
of both tha orbital and tha apla aagalar aoaaatua. Tha tlda ralaad on 
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tha aatallita by the planet will dlselpate both rotational energy and 
orbital energy. Since the rotational energy la emeh aaaller than the 
orbital energy the prooeaa will lead to reeonano* first, beoanae in 
this state, rotational energy dissipation will go to aero or to a 
snail vain* related to the changing of orbital period. Therefor*, 
starting with an eccentric orbit the dnnbbell will first settle down 
on a resonant orbit and thereafter the orbital energy will deoreas* by 
enhanced tidal dissipation on the satellite. When the eccentricity 
and senl-najor axis decreases, the systen nay shift fron a higher to a 
lower resonance. The systen will rapidly evolve if the enhanced 
dissipation Is operating. A nunerioal Integration exanple Is given 
below. The nunerioal exanple will be United to relatively snail 
eccentric ity for technical reasons. 

In principle, a dnnbbell nay be Injected on a high eccentricity 
orbit. Eventually It will lose the rotational angular nonentun and 
will stablllxe on a spin-orbital coupling reglne as does the planet 
Mercury (see Figure S-26) and as do all satellites. After being 
launched in an ecoentrlo orbit, the dnnbbell will settle In a reglne 
where It will nake. during one orbital period, as integer nuaber of 
half rotations about an axis close to the nornal to the orbital plane 
(or coinoiding with the nornal to the orbital plane if the orbit is 
equatorial or polar). Thereafter, the tidal evolution will cause the 
ecoentrioity and the senlnajor axis to decrease while the eccentricity 
is decreasing. The systen nay go through other spin orbital coupling 
resonances ending in a final stag* with the dnnbbell always aligned 
radially and noving in a circular orbit. 
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Figure 5*-26. A rough planner sketch of the orientation of 
Mercury's axis of alnlmun moment of Inertia, at different points along 
its orbit, given that the rotational period Is two-thirds of the 
orbital period and that this axis Is aligned with the Sun-planet 
vector at perihelion. 
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5.4. S Applicatioaa of System BehaviC' Orbital Trantfer 

to principle, the prooeaa Juat deeorlbed oonld be efficiently 
■echaniied in a spacecraft on a highly eccentric orbit. In 
particnlar, two examples come to mind. The first one concerns the 
transfer orbit of an OTV (see Figure S-27). One should examine if 
transfer from OEO to low Barth orbit might be accompllahed by the 
technique of tidal interaction combined perhaps with atmospheric drag. 
Starting with a transfer orbit (6800 - 42000 km) the eccentricity 
being (2a - 24400 km) e - 0.72. and a* - 24.400 (1-0.53) » 11.460. 
Without drag dissipation, this will be the radius of the circular 
final orbit after removal of energy from the orbit by tidal effects. 
The exploitation of relatively low atmospheric drag which would not 
reqnlre heavy thermal protection may drastically reduce the initial 
aemimajor axis and speed up the process of circularising at low 
altitode in the OTV parking orbit. 

The second example may be given by the capture in geostationary 

orbit of a spacecraft from a high eccentricity orbit. Naturally, if 

e 

a* ■ 42000 km. a(l-e ) has to be 42000 km also; and. since this is 
the semilatns rectum of the orbit, the perigee altitude of the 
starting orbit should be larger than 21000 km. that is, one-half of 
the minimum semilatus rectum. 

In particular, if the perigee of the orbit is 24000 km, the 
eeeentrlelty and .'‘he semimajor axis should be 

a « 96,000 km e ■ 0.75 a(l+e) » 158,000 km 

If the perigee is at 23,000 km 


a * 132,000 km e • 0.83 a(l+e) = 241,000 km 
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Figure 5-27. Orbitel transfer from GEO to LEO of a dumbbell OTV. 
Initially , the syatem may aynohronlze in a 9:1 resonance due to spin 
orbital coupling as shown. During the decay, or even from the 
beginning If the characteristics of the dumbbell and initial 
conditions are correct, the dumbbell may low order resonances - 8.5:1, 
8:1, — ——1:1. It will finally synchronize in the 1:1 spin orbital 
coupling in the final LEO orbit where the retrieval to compact 
configuration will be performed. 


i) 
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If the perigee It at 22,500 

a -148,750 e - 0.87 a(l-»-e) - 315,000 km 

It la apparent that. If the apogee should be at the Noon distance for 
tranafer from the Moon to geostationary orbit, the perigee should be 
around 22500-22600 km. 

Finally, In principle. It may be possible to Insert a dumbbell 
system In geostationary orbit by launching the complete system In a 
highly eooentrlc orbit, and either by exploitation of the Moon's 
gravity field or by a small (<200m/sec) dv, raise the perigee to the 
wanted value and bring the orbit Into the equatorial plane. 

We oomputed that without Moon assistance we still 
gain *'600*'700m/seo Av for Inserting the spacecraft In a geostationary 
orbit, ooaiparlng a classical launch from the Eastern Test Range (ETR) 
with the technique outlined above. It has to be noticed that a large 
part of the energy will be recovered during the circularization. This 
■ay supply electrical energy for the operation. Naturally one should 
be careful In considering the time required for dissipating the 
energy. It strongly depends on the design of the spring-dashpot 
system and/or of the length of the tether reeled out at perigee and 
reeled la at apogee. 

A* a technique for transferring a spacecraft from LEO to GEO with 
the aim of saving energy with respect to the energy Involved In the 
actual launch procedure from the ETR. this technique Is valid when the 
mast of the system It very large. The basic element of the transfer 
It a dissipation mechanism which absorbs energy from the orbit which 
nay be recovered as electric power. In line with what we have written 
above we want to use the enhanced effect of the tide raised by the 
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Barth oa tha apaeaoraft to daorease its aaergy. This can ba dona 


paaalvaly or aotlvaly by propar design of a aprlng-dashpot syatam and 
by proparly ooatrolllag tha tension oa tha tathar oonnaotlag tha two 
cassaa and. therefore. Its length. 


Let ns start with a spacecraft in a circular parking orbit at an 
altitude a^ (typically a^ > 220 km (equatorial) with an 

Inclination of 28** (R.T.R.) and a valoolty close to 7.8 kn/sac with 
an ineraBantal velocity 


A.v ■ V 
1 o 


-2i=-i 

a +a 

o 


(32) 


Tha spacecraft is inserted in a transfer orbit with apogee a. The 
velocity at apogee a in this Rohaan tranafer is 


V 

a 


V 

o 


21 

a +a 

o 


At apogee dlstaaoa a the circular velocity is 


V 


ca 


(33) 


(34) 


Suppose we want to raiae the perigee of the Rohaan tranafer orbit and 
change the inclination froa 28^ to 0^. We need to bring the aagnitude 
of the velocity vector to the value 


;£ iL 

i p+a 


( 35 ) 


and its orientation in the equatorial plane (see Figure 5-28). 
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Plgur* S>28. Vtotor relatlonshlpa* 
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By trlvitl oonpotatlon w« h«v« 

2a p 

A V ■ ^ 

i t(a +a) a(p-*-a) 


-4 -r^ eo$28*^ 


P*. 


1/2 


(p+i) (a +a) 
0 


(36) 


If W 9 oall a* the radius of the OBO orbit and we want to transfer the 
ipaoeoraft by tidal effect from the orbit with perigee p and apogee a 
to OEO sinoe the angular momentum is preserved a* and p have to 
satisfy the following condition 


.a - li£ 

a+p 




1+p/a 


(37) 


which meana 


'>T 


Tha relation 37) la directly dedneed from 


(38) 


a(l-e^)-a*. a(l-a)-p, a(l+o)»a 


(39) 


Tranaferrlng the apaoecraft from the orbit (a,p) to GEO implies a 
dissipation of energy 


ae - ilii. - 

2a* a-*-p 


(40) 


Hers p Is the Barth gravitational constant and M Is the mass of the 
spacecraft. If a + p *4a*) as we have shown, when we reduce the 
magnitude of iv^ Av^ or to exploit the Moon gravity assistance (a + 
p - 9a*). the amount of energy to he dissipated Is more than one-half 
of the total energy of the spacecraft in OEO. Now the total ^rbital 
energy in GEO of a spacecraft of mass N is given by Mp/2a* and, 
therefore, the order of magnitude of energy to be dissipated ranges 
between Mp/4s* and Np/2a*. A rough evaluation gives in the case of a 
+ p - 4a* 
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AE/M - 23000 Jonle/k| (41) 

Prom th«t« figures it is olssr thst it asy be diffioalt to 
lasgiae prsotiosl system wbicb ooold do the Job in s ressoasble time 
imtervsl. For exsaple. sssnalag s Nooa-ssslsted trsasfer orbit ss s 
stsrtiag orbit for sxploltstion of enhsaoed tldsl dissipstioa for 
olroalsrisstlon to OEO: the period of the first orbit will be of the 

order of 8-9 days. Siaoe the inltlsl perigee p of the orbit will be 
St roaghly 4 earth radii, at perigee passage we may depead upoa some 
amooat of eaergy disslpatloa geaerated by the tide raised at 4 earth 
radii. 


la tmQt, sappoae we start with a spacecraft of the following 
general shape (see Figare S-29) : a mala spacecraft of mass N 

ooaaected with a seooadary spacecraft with a mass of the same order of 
aagaitttde with a tether of leagth 2s. We leave for the moment the 
problem of stabilisation of the dnmbbell on a spin orbital coupling 
mode and try to estlswte the maximum work which can be done by the 
gravity gradient force. Suppose P"4E^. the grsvity field at p will be 

».••• u 


while the gravity gradient force will be 






(43) 


Since the orbit is very eccentric sad the dumbbell is very 
elongated we may assume that the angular velocity of the dnmbbell Is 
close to the orbitsl sagnlsr velocity 9 at perigee when the tide of 
the earth will keep the system oriented toward it. 
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Figure S-29. Spacecraft oonl'iguratlons for orbital pumping. 
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If n It tht mttn aotloii of the o.g, of the duabbell the orbltel 
•ngnUr Telooity tt perigee It 

, 1/2 


■ B 

P 


(1+e) 


,3/2 


(44) 


(1-e) 

whloh la the otte of the OTV-S/C It roughly 9n and la the otte 
ooatldered of t laaar ettltted orbit it of the order of 30B**31a. At 


perigee pettage the teatioa In the tether will be given by the ana of 
the gravity gradient force (43) and the differential oentrlfagal force 


Af - M 9** • M a* - Mg 

® “ ... 3 


^^(l-»e) t 
0 3 3 

a’ (1-e)’ 


(46) 


(1-e)” 

and tlnoe e It olote to 1 the tentlon will be elote to 


<«> 

Thla tenaloa It reduced to roughly 1/2 of the value (47) very rapidly 
when the gravity gradient It no longer operating and only the 
centrifugal force la keeping the tyttea In tendon. 

The work which can be done by the gravity gradient per orbit It equal 
to the elongation At of the tether at perigee multiplied by l/2r* at 
given by (47). In fact, the work to retrieve the tame height of 
tether at perigee la done agalaat the centrifugal force (due to the 
rotational motion of the dumbbell) with an angular velocity which 
remalnt almott coattant from perigee to apogee. The maximum value of 
the work dlttlpated per orbit It given by: 



or 


¥ - •• 


xAx 


( 49 ) 
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Ob th« otbar hand, tha aaouat of anatgy dlaalpated In traafarrlng 
tha apataa froa tha aeoantrio orbit to tha OBO orbit ia alaott twioa 
tha valna glvaa bj (41) that la 40000 Joala/kg. Aa avalnatlon of (49) 
glvaa 10^ Jeala/k|. not aa aaraaaoaabla valaa whan coaparad with tha 
total aaargy to ba diaslpatad. 

Thlt aaaapla ahowa that tha prooaaa nay work. In partioalar it 
nay ba aoat applieahlc to a ayataa of haga linaai diaaaaloa and if a 
ralativaly long tlaa aeala of tha prooaaa oaa ba a<*.oaptad. 

Ia ooaolaaioa, probably tha aoat iataraating practical 
applieatioaa of thaaa orbital puapiag taohniqnaa will ba for 
olroalaritlag aaall aooaatrioity orbita and for halping to circnlarlza 
an orbit prior to raaatry of an OTV in parking orbit froa Rohaan 
traaafar. Thaaa two oaaaa aaaa to daaarva particular attaution. But 
wa ahoald not oacluda tha poaaibility that farthar atudy of thaaa 
phanonana would not land to othar uaaful orbital tranafar tachniquaa. 

S.4.4 Nunarloal Sianlatloa Raaulta 

Two oaaaa wara alaaiatad naing SCYBOOC. ona with a low 
aooantrloity (a ■ 0.036 iaitlally) and aodarata altitnda (400 to 900 
ka) , tha aaeond with high aoeaatrioity (a " 0.72) and high altitnda 
(6800 ka to 42000 ka) . 

IVo aqual aaaaaa wara uaad. initially aaparatad by 100 km: tha 

SKYflOOC nodal uaad only thaaa two aaaaaa and no diaeratizatlon of tha 
tathar. In ordar to uaa eziating aoftwara for ganaratlng initial 
ooaditioaa, tha aiaalatioaa wara bagua 90^ along orbit froa porigaa 
aad ia aaatral libratioa. that ia. atraight up aad rotating at tha 
aaaa rata aa tha orbital radiaa vaotor. 
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The t«th«r oharaeterlatics «*r« detarained to that O) tha 
oattttal fraqoaacy of tha loniltudlnal tathar/aaat vlbratlona wara tha 
tana aa tha orbital pariod (l.a. raaoaaat) aad (2) tha tathar/aaaa 
ayataa waa erltloally daapad to aahaaca aaargy traaafar. 

Tha rataltt ia tha high aooantrioity orbit# although ooataining 
faataraa aipaotad froa tha thaory. wara coafuaiag aad dua to tiaa aad 
fuadiag ooaatraiata. wa wara anabla to fally iavattigata tha oaatat 
for thaaa iaooaaiataat raaalta. Mott likaly, tha problea liat ia tha 
way this run waa tat up tinea tha agraaaant batwaan tha thaory and tha 
aodal waa axoallant ia tha low acoantrieity oata. 

Tha ratultt of tha low aoeantrioity ttudy. whioh ran for about 22 
orbitt, ara oontiataat with tha thaory ontlinad abova. Prom tha 
SKYHOOK output, a varitty of paraaatart of tha 6aatar'*of-«att orbit, 
at wall at tha tathar laagth and llbration angla of tha tyatam. wara 
eonputad and printad. 

Tha aaaantrioity aad taai*aaJor aria thowad tba azpaetad 
bahavior: a pariodio (at tha orbital pariod) eonpoaaat tuparlmpotail 

on a taaular daaraata. Tha aeoaatriaity it thowa ia Pigura S-30, 
whtra both tha pariodio ad taaular tarma ara olaar. 

Por a aiora quaatitativa diteuttioa. aoatidar tba ratultt for ona 
partioular orbit. Apply aqaatloa 11) abova, Aa/orbit ■ -9 Az.z/a , 
and aompara thit to tha obaarvad ehanga ia aoeantrioity. For orbit 5. 
tha tathar langth (taparation of nataata) it thowa in Figure S-31a, 
tha miaiaum aad maaiaua ara 99. 8d aad 108.22 km. Tha aeeaatrioity at 
aaaeaativa paaka sra 0.035142 aad 0.034994, glviag Aa ■ -1.48 x 10~^. 
To apply 11). raoall that t it tha tathar langth whila At it tha 

full aaonat of tha iaoraata or daaraata. That wa uaa tha avaraga 2z ■ 
(89. 8d ^ 108.22)/2 or z • 49.52 km. aad Az - 108.22 - 89.86 - 18.36 


■ r 
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Fleure S~30. Eoocntrlcity of the center-of^Bass orbit vs. the 
orbital phase (true anosaly/ln) . A critically dasped dumbbell 
oonflfuration. with ooclllatory frequenoy tuned to the orbit, la 
allowed to evolve, converting orbital energy to internal (e.g. 
thermal) energy, oiroularlsing the orbit. The nunerioal aiaulatlon 
wee made with SKYHOOK. 
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Figure 5-31 (a). Tether length ve. orbital phase for one 
particular orbit of the study shown in Figure 3*30. (b) The 

derivative of the curve shown in (a), together with the gravity 
gradient forcing function. These are nearly in phase and their 
product (the work done by the orbit on the dumbbell) Is nearly 
maximal, (o) A plot similar to (b) for a case with zero time lag in 
tether response. The work is shown as well, and obviously Integrates 
to zero over an orbit. 








FORCING FUNCTION (Gravity Gradient) 


d (Tether Length) 
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km. With « aeml-major tzlt of 7028. 7S km (the mean of the maximum and 
minimum for the orbit), this give* an expected de/orbit of 1.66 x 


Thla agreement la qnite good. Indeed anrprlalngly good 
oonildarini that the predicted value is for a maximally efficient 
tyatem in which all the work it done on or by the ayatem at momenta of 
maximal or minimal tenaion. However, in tome aente the critical 
damping oantet the dumbbell to behave "well": the work done over an 
orbit la more properly / . (Force) dt. Figure 5-31(b) ahowa 

the gravity gradient forcing function (arbitrarily normalized and with 
the d.o. component, which integratea to zero, removed) together with 
the obaerved d(2z)/dt. Theae two are nearly in phaae, ao their 
product ia (almoat) alwayt poaitive and the Integral, hence the work 
done by the orbit, ia (nearly) optimal with reapect to the tether* a 
reaponaa lag. A aero-lag ideallxed oaae ia ahown in Fignre S-31(o); 
the integral over an orbit ia aeen to be zero. 

We can compute, to firat order, thetiiru.* required to circularize 

the orbit aa followa. Since Ae/orbit it tmall, we may write 11) aa 
2 

de/dN • -97 Az/a . where N ia the orbit number. If we aatnme that, 
even allowing for the phaae lag. Az it proportional to the gravity 
gradient forcing function (g.g.f.), then (to firat order) Az 

• [max-mini g.g.f. « [max-mini (1/r^) « [max-mini 1/ (l-t-e cos9)^ 

* [max-mini (l-3ecot6)« e. Hence Az ■ Az (e/e ) where Az ,e and 

0 0 0 0 

2 

z^ are the values at N • 0, aay. Then de/dN « [-9z Az /(a e )1e » 

0 0 0 0 

((de/dN) (1/e )la, and e(N) ■ e exp(-N/t) where t * e /(de/dN) is 

a "time" scale for decay. In the low eccentricity case, from the 
observed Ae/orbit, we have t ■ 3.5 z 10~^/1.S x lO”^ ■ 230 orbits, or 
about IS days. To oircnlarize an orbit from, say, S% to 1% 
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eoocatrioity would that take about log (O.OS/O.ODt ~ 9.S dayt. a 
praotleal and utaful length of time for such an operation. 

The noaerloal elanlation reaulta therefore oonfim the theory, at 
leaat in the low eooentrioity case, and the ntefulneta of orbital 
puaping for orbit olroularization is deaonatrated in a practioal 

aense. 
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